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The  experimental  results  described  in  this  report  are  concerned  with  a 
nominally  two-dimensional  separating  turbulent  boundary  layer  for  an  airfoil- 
type  flow  in  which  the  flow  was  accelerated  and  then  decelerated  until  sepa- 
ration„  Upstream  of  separation  single  and  cross-wire  hot-wire  anemometer 
measurements  are  also  presented.  Measurements  obtained  in  the  separated  zone 
with  a  directionally-sensitive  laser  anemometer  system  are  presented, for  U, 

V,  u2,  v2,  -uv,  u3,  u4,  v3,  v4,  the  fraction  of  time  that  the  flow  moves  down-  j 
stream,  the  fraction  of  time  that  the  flow  moves  away  from  the  wall,  and  u 
spectra. 

In  addition  to  confirming  the  earlier  conclusions  of  Simpson  et  al .  (1977) 
regarding  a  separating  airfoil-type  turbulent  boundary  layer,  much  new  infor¬ 
mation  about  the  separated  region  has  been  gathered.  (1)  The  backflow  mean 
velocity  profile  scales  on  the  maximum  negative  mean  velocity  Ur  and  its  dis¬ 
tance  from  the  wall  N.  A  U+  vs.  y+  law-of-the-wall  velocity  profile  is  not 
consistent  with  this  result.  (2)  The  turbulent  velocities  are  comparable  with 
the  mean  velocity  in  the  backflow,  although  low  turbulent  shearing  stresses 
are  present.  (3)  Mixing  length  and  eddy  viscosity  models  are  physically 
meaningless  in  the  backflow.  (4)  Negligible  turbulence  energy  production  occurs 
in  the  backflow. 

V' 

These  and  other  results  lead  to  significant  conclusions  about  the  nature 
of  the  separated  flow  When  the  thickness  of  the  backflow  region  is  small  as 
compared  with  the  shear  layer  thickness.  The  backflow  is  controlled  by  the 
large-scale  outer  region  flow.  The  small  mean  backflow  does  not  come  from  far 
downstream,  but  appears  to  be  supplied  intermittently  by  large-scale  structures 
as  they  pass  through  the  separated  flow.  Downstream  of  fully-developed  sepa¬ 
ration,  the  mean  backflow  appears  to  be  divided  into  three  layers:  a  viscous 
layer  nearest  the  wall  that  is  dominated  by  the  turbulent  flow  unsteadiness  but 
with  little  Reynolds  shearing  stress  effects;  a  rather  flat  intermediate  layer 
that  seems  to  act  as  an  overlap  region  between  the  viscous  wall  and  outer  re¬ 
gions;  and  the  outer  backflow  region  that  is  really  part  of  the  large-scaled 
outer  region  flow.  .The  Reynolds  shearing  stress  must  be  modeled  by  relating 
it  to  the  turbulence!  structure  and  not  to  local  mean  velocity  gradients.  The 
mean  velocities  in  rhe  backflow  are  the  results  of  time-averaging  the  large 
turbulent  fluctuations  and  are  not  related  to  the  source  of  the  turbulence. 
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ABSTRACT 

The  problem  of  turbulent  boundary  layer  separation  due  to  an 
adverse  pressure  gradient  is  an  old  but  still  important  problem  in 
many  fluid  flow  devices.  Until  recent  years  little  quantitative  ex¬ 
perimental  information  was  available  on  the  flow  structure  downstream 
of  separation  because  of  the  lack  of  proper  instrumentation.  The 
directionally-sensitive  laser  anemometer  now  provides  the  ability  to 
accurately  measure  the  instantaneous  flow  direction  and  magnitude. 

The  experimental  results  described  in  this  report  are  concerned 
with  a  nominally  two-dimensional  separating  turbulent  boundary  layer 
for  an  airfoil-type  flow  in  which  the  flow  was  accelerated  and  then 
decelerated  until  separation.  Upstream  of  separation  single  and  cross¬ 
wire  hot-wire  anemometer  measurements  are  also  presented.  Measurements 

obtained  in  the  separated  zone  with  a  directionally-sensitive  laser 
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anemometer  system  are  presented  for  U,  V,  u  ,  v  ,  -uv,  u  ,  u  ,  v  ,  v  , 

the  fraction  of  time  that  the  flow  moves  downstream,  the  fraction  of  time 

that  the  flow  moves  away  from  the  wall,  and  u  spectra. 
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In  addition  to  confirming  the  earlier  conclusions  of  Simpson 
et  al.  (1977)  regarding  a  separating  airfoil-type  turbulent  boundary 
layer,  much  new  information  about  the  separated  region  has  been  gath¬ 
ered.  (1)  The  backflow  mean  velocity  profile  scales  on  the  maximum 
negative  mean  velocity  UN  and  its  distance  from  the  wall  N.  A  U+  vs.  y+ 
law- of- the- wall  velocity  profile  is  not  consistent  with  this  result. 

(2)  The  turbulent  velocities  are  comparable  with  the  mean  velocity  in 
the  backflow,  although  low  turbulent  shearing  stresses  are  present. 

(3)  Mixing  length  and  eddy  viscosity  models  are  physically  meaningless 
in  the  backflow.  (4)  Negligible  turbulence  energy  production  occurs  in 
the  backflow. 

These  and  other  results  lead  to  significant  conclusions  about  the 
nature  of  the  separated  flow  when  the  thickness  of  the  backflow  region 
is  small  as  compared  with  the  shear  layer  thickness.  The  backflow  is 
controlled  by  the  large-scale  outer  region  flow.  The  small  mean  back- 
flow  does  not  come  from  far  downstream,  but  appears  to  be  supplied  inter¬ 
mittently  by  large-scale  structures  as  they  pass  through  the  separated 
flow.  Downstream  of  fully-developed  separation,  the  mean  backflow  appears 
to  be  divided  into  three  layers:  a  viscous  layer  nearest  the  wall  that 
is  dominated  by  the  turbulent  flow  unsteadiness  but  with  little  Reynolds 
shearing  stress  effects;  a  rather  flat  intermediate  layer  that  seems  to 
act  as  an  overlap  region  between  the  viscous  wall  and  outer  regions;  and 
the  outer  backflow  region  that  is  really  part  of  the  large-scaled  outer 
region  flow.  The  Reynolds  shearing  stress  must  be  modeled  by  relating  it 
to  the  turbulence  structure  and  not  to  local  mean  velocity  gradients. 

The  mean  velocities  in  the  backflow  are  the  results  of  time-averaging  the 
large  turbulent  fluctuations  and  are  not  related  to  the  source  of  the  tur¬ 
bulence. 
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I.  INTRODUCTION 

The  problem  of  turbulent  boundary  layer  separation  due  to  an  adverse 
pressure  gradient  is  an  important  factor  in  the  design  of  many  devices  such 
as  jet  engines,  rocket  nozzles,  airfoils  and  helicopter  blades,  and  the  design 
of  fluidic  logic  systems.  Until  the  last  five  years  little  quantitative  experi¬ 
mental  information  was  available  on  the  flow  structure  downstream  of  separation 
because  of  the  lack  of  proper  instrumentation. 

In  1974  after  several  years  of  development,  a  one  velocity  component 
directionally-sensitive  laser  anemometer  system  was  used  to  reveal  some  new 
features  of  a  separating  turbulent  boundary  layer  (Simpson  et  al . ,  1974). 

The  directional  sensitivity  of  the  laser  anemometer  system  was  necessary  since 
the  magnitude  and  direction  of  the  flow  must  be  known  when  the  flow  moves  in 
different  directions  at  different  instants  in  time.  In  addition  to  much  turbu¬ 
lence  structure  information,  it  was  determi ned: ( 1 )  that  the  law-of-the-wall 
velocity  profile  is  apparently  valid  up  to  the  beginning  of  intermittent  separa¬ 
tion:  (2)  that  the  location  of  the  beginning  of  intermittent  separation  or  the 
upstreammost  location  where  separation  occurs  intermittently  is  located  close 
to  where  the  freestream  pressure  gradient  begins  to  rapidly  decrease;  (3)  that 
the  normal  stress  terms  of  the  momentum  and  turbulent  kinetic  energy  equations 
are  important  near  separation;  and  (4)  that  the  separated  flowfield  shows  some 
similarity  of  the  streamwise  velocity  U,  of  the  velocity  fluctuation  u',  and  of 
the  fraction  of  time  that  the  flow  moves  downstream  (Simpson  et  al . ,  1977). 

Based  upon  these  results,  modifications (Simpson  and  Collins,  1978;  Collins 
and  Simpson,  1978)  to  the  Bradshaw  et  al .  (1967)  boundary  layer  prediction 
method  were  made  with  significant  improvements.  However,  this  prediction  effort 

pointed  to  the  need  to  understand  the  relationship  between  the  pressure  gradient 


1 


relaxation  and  the  intermittent  separation  region  structure.  A  number  of 
other  workers  have  tried  to  predict  this  type  flow,  but  with  questionable 
assumptions  about  the  turbulence  structure  near  the  wall.  In  nearly  all 
efforts,  the  workers  have  simply  extended  the  velocity  and  turbulence  profile 
correlations  that  apply  to  attached  flows  to  the  backflow  region.  Even  though 
turbulent  fluctuations  near  the  wall  in  the  backflow  region  are  as  large  as  or 
larger  than  mean  velocities,  these  predictors  use  a  turbulence  model  that  is 
tied  to  the  mean  velocity  gradient.  Even  with  adjustment  of  turbulence  model 
"constants"  to  fit  one  feature  or  another,  these  models  do  not  predict  simul¬ 
taneously  the  backflow  velocity  profile,  the  steamwise  pressure  distribution, 
and  the  fact  that  length  scales  increase  along  the  flow.  Clearly  then,  a 
limiting  factor  for  further  improvement  of  the  prediction  of  separated  flows 
is  the  lack  of  fundamental  experimental  velocity  and  turbulence  structure  infor¬ 
mation  with  which  to  develop  adequate  models,  especially  for  the  backflow  region. 
Such  data  are  presented  here. 

The  experimental  results  described  in  this  report  are  concerned  with  a 

nominally  two-dimensional  separating  turbulent  boundary  layer  for  an  airfoil- 

type  flow  in  which  the  flow  was  accelerated  and  then  decelerated  until  separation. 

Upstream  of  separation  single  and  cross-wire  hot-wire  anemometer  measurement 

results  are  presented.  Measurements  obtained  in  the  separated  zone  with  a  direc- 

~2  ~~2  _ 

tionally-sensitive  laser  anemometer  system  are  presented  for  U,  V,  u  ,  v  ,  -uv, 
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u  »  v  ,  u4,  v  ,  the  fraction  of  time  that  the  flow  moves  downstream  yu>  fraction 
of  time  that  the  flow  moves  away  from  the  wall  y  ,  and  u  spectra.  The  implications 
of  these  results  to  flow  models  are  discussed. 

II  EXPERIMENTAL  EQUIPMENT 


II. 1  Basic  Wind  Tunnel 


, .  -..,  ....... 


The  mainstream  flow  of  the  blown  open-circuit  wind  tunnel  is  intro¬ 
duced  into  the  test  section  after  first  passing  through  a  filter,  blower,  a 
fixed-setting  damper,  a  plenum,  a  section  of  honeycomb  to  remove  the  mean 
swirl  of  the  flow,  seven  screens  to  remove  much  of  the  turbulence  intensity, 
and  finally  through  a  two-dimensional  4:1  contraction  ratio  nozzle  to  further 
reduce  the  longitudinal  turbulence  intensity  while  accelerating  the  flow  to 
test  speed.  These  same  components  were  in  an  earlier  version  of  this  wind 
tunnel  with  a  shorter  test  section  that  was  used  in  previous  research  (Simpson 
et  al . ,  1977;  Simpson  and  Wallace,  1975;  Simpson  and  Shackleton,  1977). 

Fig.  1  is  a  side  view  schematic  of  the  25  feet  long,  3  feet  wide  test 
section  of  the  wind  tunnel.  The  upper  wall  is  adjustable  such  that  the  free- 
stream  velocity  or  pressure  gradient  can  be  adjusted.  The  side  walls  are  made 
of  plexiglass.  The  test  wall  is  constructed  from  3/4  inches  thick  fin-form 
plywood,  reinforced  every  11  inches  with  3  x  1%  x  %  inches  cross  section  steel 
channel . 

The  active  boundary  layer  control  system,  which  is  described  in  section 
I I. 2  below,  is  used  to  eliminate  preferential  separation  of  the  curved  top  wall 
boundary  layer.  Highly  two-dimensional  wall  jets  of  high  velocity  air  are 
introduced  at  the  beginning  of  each  of  the  eight  feet  long  sections.  At  the 
latter  two  streamwise  locations  the  oncoming  boundary  layer  is  partially  removed 
by  a  highly  two-dimensional  suction  system.  In  order  to  accorrenodate  the  increased 
energy  dissipation  in  the  wind  tunnel  laboratory  due  to  the  boundary  layer  control 
system,  a  new  3  ton  air  conditioner  was  added. 

The  inviscid  core  flow  is  uniform  within  0.05%  in  the  spanwise  direction 
and  within  1%  in  the  vertical  direction  with  a  turbulence  intensity  of  0.1%  at 
60  fps.  The  test  wall  boundary  layer  is  tripped  by  the  blunt  leading  edge  of  the 
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Figure  1.  Sideview  schematic  of  the  test  section.  Major  divisions  on  scales: 
10  inches.  Note  baffle  plate  upstream  of  blunt  leading  edge  on 
bottom  test  wall  and  side  and  upper  wall  jet  boundary  layer  controls. 
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Figure  2.  Cross  sectional  view  of  wall  jet-wall  suction  assembly. 
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plywood  floor,  the  height  of  the  step  from  the  wind  tunnel  contraction  to  the 
test  wall  being  h  inches.  Smoke  can  be  introduced  uniformly  into  the  boundary 
layer  just  upstream  of  this  trip  for  use  with  the  laser  doppler  anemometer. 

II .2  Boundary  Layer  Control  System 

An  active  boundary  layer  control  system  was  installed  on  the  nontest  walls 
of  the  test  section  to  inhibit  undesirable  flow  three-dimensionality  and  to  pre¬ 
vent  separation.  Because  the  static  pressure  in  the  test  section  is  time  varying 
in  a  series  of  unsteady  experiments  at  the  same  mean  conditions  as  in  the  case 
reported  here  (Simpson  et  al . ,  1980),  no  passive  boundary  layer  control  can  be 
used  that  depends  on  a  steady  test  section  pressure  higher  than  the  pressure  out¬ 
side  the  tunnel.  In  previous  steady  freestream  separated  flow  experiments  in  the 
old  tunnel  test  section  (Simpson  et  al . ,  1977),  a  perforated  plate  was  located  at 
the  test  section  exit  to  produce  static  pressures  in  the  test  section  that  were 
above  atmospheric  pressure.  In  that  case,  the  upper  wall  boundary  layer  was  bled 
from  the  test  section  just  upstream  of  the  separation  zone. 

In  the  current  case,  the  active  boundary  layer  control  system  removes  low 
momentum  fluid  by  sucking  off  the  boundary  layer  and  supplies  high  momentum  fluid 
through  trangential  wall  jets.  Its  performance  is  less  influenced  by  the  fluctu¬ 
ating  test  section  static  pressure  in  the  unsteady  experiments  than  that  of  the 
previous  passive  system.  Fig.  2  is  a  cross-sectional  view  of  the  wall  suction  and 
wall  jet  units  located  at  100  inches  and  200  inches  on  the  non-test  walls.  Only 
the  wall  jet  portion  of  this  unit  is  installed  at  the  test  section  entrance.  Only 
the  essential  features  of  this  system  are  summarized  here;  other  details  are  con¬ 
tained  in  Simpson  et  al .  (1980)  and  in  an  unpublished  report  by  Bowles  (1977). 

A  fan  supplies  the  wall  jets  with  air  sucked  through  the  suction  slots  and 
with  additional  makeup  air.  Sheet  metal  ducting,  flow  dampers,  resistance  baffles, 
and  manifolds  that  contain  the  wall  suction  and  wall  jet  units  are  used  to  direct 
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and  control  the  air  distribution.  All  of  the  wall  suction  and  wall  jet 
units  have  identical  cross  sections.  As  much  care  as  possible  was  taken 
to  make  these  units  geometrically  and  aerodynamically  two-dimensional.  As 
shown  in  Fig.  2  the  wall  jet  portion  is  a  6:1  area  ratio  nozzle  that  accelerates 
the  fluid  before  it  is  injected  along  the  glass  tunnel  sides  walls  or  plexiglas 
top.  The  suction  portion  is  a  1:3.6  area  ratio  diffuser  that  decelerates  the 
removed  flow.  An  aluminum  divider  plate  separates  the  wall  jet  and  suction 

flows,  forming  one  wall  of  the  jet  nozzle  and  one  wall  of  the  diffuser.  The 

specially  extruded  aluminum  deflector  directs  the  jet  flow  parallel  to  the 
tunnel  wall  and  scoops  the  suction  flow  from  the  upstream  tunnel  flow.  Shims 
were  placed  between  the  aluminum  divider  plate  and  the  extruded  aluminum  deflector 
to  make  the  jet  flow  exit  gap  uniform  within  0.0016  inches  along  a  given  unit. 

The  gap  was  nominally  0.25  inches,  but  was  slightly  different  for  each  unit.  The 

suction  side  flow  entrance  gap  was  less  uniform,  being  within  0.01  inches  of  the 

nominal  5/16  inches. 

As  also  shown  in  Fig.  2,  32  mesh  stainless  steel  screen  and  1/8  inches  cell 
hexagonal  aluminum  honeycomb  are  located  at  the  nozzle  entrance  to  evenly  distri¬ 
bute  and  straighten  the  flow  from  the  jet  manifold.  After  initial  system  testing, 
additional  screen  was  mounted  on  top  of  the  honeycomb  opposite  the  manifold  duct 
entrance  from  the  supply  duct.  This  eliminated  preferential  flow  due  to  infringe¬ 
ment  of  the  incoming  flow.  In  one  case  a  flow  deflector  was  required  at  the 
manifold  entrance  to  further  distribute  the  flow.  32  mesh  screen  was  also  placed 
over  the  diffuser  exit  to  distribute  the  suction  flow  more  evenly  over  the  flow 
cross-sectional  area.  A  small  gap  between  the  screen  and  each  endplate  was  used 
to  induce  a  greater  flow  along  the  endplates.  This  greater  flow  benefits  the 
momentum  deficient  wind  tunnel  corner  flows. 
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Two  pressure  taps  were  located  in  the  nozzle  and  in  the  diffuser. 

After  calibration,  the  measured  pressure  difference  between  these  taps 
allowed  the  nozzle  and  the  diffuser  to  be  used  as  jet  flow  and  suction 
flow  meters,  respectively.  Excellent  linear  calibrations  were  found 
between  measured  dynamic  pressures  and  the  respective  differential 
pressures. 

The  average  dynamic  pressure  of  the  jet  exit  flow  was  measured  along 
the  length  of  each  unit  with  a  0.25  inches  dia.  impact  probe.  The  standard 
deviation  of  the  dynamic  pressure  variation  was  less  than  2.5%  along  each  of 
the  jet  units.  The  dynamic  pressure  in  the  3/4  inches  nearest  the  end  of 
each  unit  was  about  2/3  of  that  for  the  midsection.  The  jets  at  the  test 
section  entrance  were  operated  at  an  average  velocity  of  90  fps;  at  the  100 
inches  location  the  wall  jet  velocity  was  120  fps  for  the  upper  wall  and  72 
fps  for  the  side  walls;  at  the  200  inches  location  the  upper  wall  and  side 
wall  jet  velocities  were  75  fps  and  57  fps. 

Fig.  3  shows  the  mean  velocity  and  streamwise  turbulence  intensity 
profiles  in  the  midplane  along  the  second  streamwise  upper  wall  jet  that  were 
obtained  from  a  hot-wire  anemometer.  Note  that  the  velocity  profile  is  asymmetric 
with  more  high  velocity  flow  near  the  freestream  side  of  the  jet.  Saripalli  and 
Simpson  (1979)  have  shown  that  such  an  asymmetric  jet  is  more  effective  in  pre¬ 
venting  boundary  layer  separation  than  a  uniform  jet  with  a  greater  momentum  flow 
rate.  This  is  due  to  the  fact  that  less  of  the  asymmetric  jet  momentum  is  wasted 
on  increased  wall  drag  while  greater  mixing  occurs  with  the  outer  region  flow. 

The  variation  of  the  dynamic  pressure  of  the  suction  flow  was  measured 
along  the  length  of  each  unit.  The  difference  between  the  static  pressure  inside 
the  diffuser  at  a  particular  location  and  atmospheric  pressure  had  a  standard 
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deviation  of  no  more  than  2%  along  each  suction  unit.  Most  of  the 
departure  from  uniformity  occurred  near  the  ends  where,  fortunately, 
about  20%  greater  suction  dynamic  pressure  occurred.  Thus  a  greater 
amount  of  momentum  deficient  wind  tunnel  corner  fluid  could  be  removed, 
partly  overcoming  the  effect  of  the  dynamic  pressure  deficiency  at  the 
ends  of  the  wall  jet  sections.  The  average  suction  velocity  at  the  100 
inches  location  was  45  fps  while  for  the  200  inches  location  it  was  75 
fps  for  the  upper  wall  and  56  fps  for  the  side  walls. 

The  hot-wire  anemometer  mean  velocity  and  streamwise  turbulence 
intensity  profiles  in  the  midplane  along  the  first  streamwise  upper  wall 
suction  unit  were  almost  uniform.  This  indicates  that  some  of  the  upstream 
flow  is  deflected  toward  the  wall  by  the  suction  since  the  boundary  layer 
velocity  profile  is  not  uniform.  Immediately  downstream  the  asymmetric 
velocity  profile  jet  momentum  is  then  rapidly  mixed  with  the  remaining  up¬ 
stream  boundary  layer  flow. 

It  should  be  noted  that  the  flows  in  this  control  system  are  relatively 
insensitive  to  the  +1  cm  of  water  static  pressure  oscillation  in  the  test 
section  in  the  unsteady  experiments.  The  large  volume  of  the  control  system 
and  the  12  inches  of  water  static  pressure  loss  in  its  components  act  as  a 
large  low-pass-frequency  filter.  Dynamic  pressure  oscillations  of  the  wall 
jet  flow  were  of  the  order  +  0.016  inches  of  water. 

1 1. 3  Hot-wire  Anemometers 

Miller-type  (1976)  integrated  circuit  hot-wire  anemometers  and  linearizers, 
as  modified  by  Simpson  et  al .  (1979)  were  constructed  and  used.  A  TSI  Model  1050 
anemometer  was  used  with  the  surface  hot-wire  element  that  is  described  in  section 

11. 4  below.  The  frequency  response  was  flat  up  to  7.5  kHz  for  an  overheat  ratio 
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of  0.7.  this  moderately  high  overheat  ratio  was  used  since  Wood  (1975) 
showed  that  the  range  of  flat  frequency  response  is  improved  with  a  higher 
overheat  ratio. 

Standard  TSI  model  1274-TI.4  normal  wire  and  model  1248-TI.5  cross¬ 
wire  probes  were  used  for  boundary  layer  measurements.  The  closest  to  the 
wall  that  these  probes  could  safely  make  measurements  was  about  0.002  inches 
and  0.035  inches,  respectively.  The  sensing  elements  are  0.00015  inches 
diameter,  0.050  inches  length  platinum-plated  tungsten  wires. 

The  traversing  mechanism  used  for  the  boundary  layer  velocity  measure¬ 
ments  was  mounted  on  the  supporting  frame  for  the  upper  wall  and  provided  for 
precise  positioning  of  the  probe  sensors  as  described  by  Strickland  and  Simpson 
(1973).  A  cathetometer  was  used  to  accurately  locate  the  probe  sensor  from  the 
wall  within  an  uncertainty  of  about  +  0.002  inches.  The  detailed  streamwise 
free-stream  velocity  distributions  were  obtained  using  a  the  Model  1274-TI.5 
probe  mounted  on  the  toy  racing  car  shown  in  Fig.  3  of  Simpson  and  Wallace  (1975). 
The  car  was  easily  positioned  along  the  flow  by  fishing  line. 

Calibrations  were  made  in  a  TSI  Model  1127  calibrator.  There  was  no 
detectable  drift  of  the  anemometer;  the  function-module  type  linearizers  had 
a  small  mount  of  DC  drift.  Each  linearized  calibration  had  a  low  level  of 
dispersion  from  a  straight  line,  with  a  product  moment  correlation  coefficient 
(Bragg,  1974)  in  excess  of  0.09999.  The  slope  of  each  calibration  varied  no 
more  than  about  4%  over  the  life  of  a  given  probe. 

A  standard  TSI  model  1015C  correlator  was  used  to  obtain  sum  and  difference 

values  for  u  and  v  from  cross-wire  signals.  Electronic  multipliers  (Analog 

2  3 

Devices  AD533JH)  were  used  to  produce  the  turbulence  quantities  uv,  u  v,  and  v  . 
These  were  trimmed  to  within  +  1%  full  scale  nonlinearity  error.  True  inte- 


grating  voltmeters,  consisting  of  a  voltage-controlled  oscillator 
(Tektronix  F6501  Function  Generator)  and  a  digital  counter  (Tektronix 
DC503  Universal  Counter),  were  used  to  obtain  true  time-averaged  results. 

I I. 4  Surface  Hot-wire  Skin  Friction  Gage 

Because  a  single  universal  calibration  is  valid  in  both  laminar  and 
turbulent  flow  and  is  insensitive  to  pressure  gradients  (Murthy  and  Rose,  1978; 
Higuchi  and  Peake,  1978),  a  Rubesin  et  al .  (1975)  type  surface  hot-wire  skin 
friction  gage  was  constructed  and  used.  The  basic  advantages  of  this  type  gage 
are  that  the  surface-heating-element  dimension  in  the  streamwise  direction  is 
very  small  and  that  conduction  losses  to  a  very  low  thermal  conductivity  substrate 
are  minimized. 

A  0.001  inches  diameter  platinum  -  10%  rhodium  wire  was  mounted  between 
0.052  inches  diameter  nickel  electrodes  located  0.4  inches  apart  whose  ends  were 
flush  with  the  flat  polystyrene  surface.  Conduction  losses  to  the  electrodes  are 
small  since  the  wire  length-to-diameter  ratio  of  400  is  large.  Several  drops  of 
ethyl  acetate  were  used  to  dissolve  the  polystyrene  in  the  vicinity  of  the  wire 
and  imbed  it  in  the  surface.  The  ends  of  the  wire  were  then  soldered  to  the 
electrodes.  A  0.00015  inches  diameter  wire  was  tried  but  was  too  fragile  for  use 
with  simple  construction  techniques.  The  polystyrene  was  mounted  on  a  thin 
portable  plexiglas  plate.  The  resulting  surface  was  sanded  and  polished  flat  and 
smooth  before  the  wire  was  mounted.  This  plate  allows  a  single  element  to  be 
moved  to  variaous  measurement  locations  with  a  minimum  of  flow  disturbance.  The 
element  is  sufficiently  downstream  of  the  end  of  the  small  ramp  and  sufficiently 
upstream  of  the  trailing  edge  to  avoid  sensing  local  disturbances  generated  by 
the  plate.  A  0.001  inches  diameter  platinum  and  platinum-10%  rhodium  thermo¬ 
couple  was  mounted  3/32  inches  downstream  of  the  hot-wire  element. 
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Rubesin  et  al .  found  that  overheat  temperatures  of  at  least  80°F  were 

needed  to  make  the  heat  loss  from  a  wire  proportional  to  its  temperature  rise, 

2 

or  E  /RAT  a  constant.  Peake  and  Higuchi  found  that  overheats  greater  than 
176°F  caused  the  wire  to  melt  the  substrate  and  separate  from  the  surface. 

Here  the  cold  resistance  at  77°F  was  3.70  ft  and  0.5  ft  overheat  resistance  was 
used,  so  with  a  temperature  coefficient  of  resistivity  of  0.89  x  10"^  °F 
then  AT  was  152°F.  The  wire  was  not  observed  to  separate  from  the  surface. 

A  simple  stainless  steel  cone  with  0.5°  angle  between  the  cone  and  the 

plate  surface  was  constructed  for  calibration  of  this  gage.  A  brass  housing 

held  the  cone  in  place  on  the  plate.  The  hot-wire  was  aligned  with  a  radial 

line  from  the  cone  apex.  The  velocity  gradient  at  the  plate  surface  was  inde¬ 

pendent  of  the  radial  position  since  the  cone  surface  velocity,  wr,  and  the 
spacing  between  the  cone  and  the  plate,  r  tan  (0.5°),  each  vary  linearly  with 
the  radius.  Because  the  maximum  surface  velocity  gradient  of  interest  was 
about  9.6  x  10  sec,  a  high-speed  grinder  motor  (Sears  and  Roebuck,  Model 
315.17440,  26000  RPM)  and  a  Variac  power  control  were  used  to  produce  600  rpm 
<  f  <  8000  rpm.  A  vinyl  tubing  flexible  connector  was  used  between  the  cone 
shaft  and  the  grinder  to  minimize  misalignment.  The  angular  speed  f  was  measured 
by  reflecting  a  light  from  the  hexagonal  grinder  chuck  nut  into  a  photomultipler 
tube  and  counting  the  signal  pulse  rate  fp  on  a  digital  counter;  thus  f  =  f p/6 . 
Heating  of  the  calibrator  flow  occurred  above  8000  rpm  due  to  substantial  frictional 
heating  in  the  steel -brass  bearing.  Since  the  air  temperature  was  measured  with 
the  thermocouple,  corrections  could  be  made.  After  calibration,  a  Mi ller- type 
exponential  electronic  linearizer  was  used  to  linearize  the  bridge  output  voltage. 

II .5  Laser  Anemometer  and  Signal  Processing 

The  laser  anemometer  used  in  these  experiments  and  shown  in  Fig.  4  is 
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Figure  4.  Photograph  of  current  two-velocity  component  di¬ 
rectional  ly-sensitive  fringe-type  laser  anemometer. 
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described  in  some  detail  by  Simpson  and  Chew  (1979).  In  essence  this  is 
a  two-velocity-component  (l),  V)  directionally-sensitive  fringe  type  system 
that  has  been  used  in  earlier  work  (Simpson  et  al . ,  1977).  The  unshifted 
and  25  MHz  Bragg-cell  shifted  beams  lie  in  an  almost  horizontal  plane  and 
measure  the  streamwise  velocity  with  vertical  fringes.  The  unshifted  and 
15  MHz  Bragg-cell  shifted  beams  lie  in  a  vertical  plane  and  measure  (v  cos  4.4° 

+  v  si n  4.4°)  with  almost  horizontal  fringes.  The  25  MHz  and  15  MHz  beams  form 
a  third  fringe  pattern  that  measures  ( U  -  T  cos  4.4°  -  W  sin  4.4°)/y£~around 
10  MHz.  Since  u^  and  (v  cos  4.4°  +  w  sin  4.40)2  were  measured  independently 
and  uw  was  presumed  very  small,  the  Reynolds  shearing  stress  -uv  resulted  from 
this  measurement.  Signal  processing  was  by  fast-sweep-rate  sampling  spectrum 
analysis,  as  described  by  Simpson  and  Barr  (1975). 

The  1  micron  dioctal  phthalate  particles  follow  the  highly  turbulent 
oscillations  found  in  separated  regions  (Simpson  and  Chew,  1979).  It  should  be 
noted  that  it  is  impossible  to  seed  a  highly  turbulent  flow  in  any  prescribed 
manner.  Highly  turbulent  flows  are  characterized  by  intense  mixing  with  the 
flow.  In  this  case  there  is  also  significant  entrainment  of  freestream  fluid 
into  the  turbulent  motions.  This  would  progressively  dilute  the  particle 
concentration  if  only  the  shear  flow  has  been  seeded.  Instead  of  needless 
worry  over  prescribed  particle  concentration,  we  have  been  concerned  with  proper 
averaging  of  available  signals  as  described  below,  with  enough  particles  to 
provide  a  high  data  rate,  and  with  sufficiently  small  particles  to  accurately 
follow  the  flow.  In  fact,  without  any  seeding  we  were  able  to  obtain  signals 
from  ambient  dust.  However,  we  used  minimal  seeding  to  produce  a  signal  data 
rate  of  about  400  per  second. 

Since  the  particle  number  density  in  a  highly  turbulent  flow  cannot  be 


made  uniform,  the  time  between  the  passage  of  successive  signal  generating 
particles  will  be  unequal.  This  effect  alone  presents  no  particular  signal 
processing  problem  if  the  time  intervals  between  successive  signal  bursts  are 
small  compared  to  l/fmax»  the  time  period  of  the  highest  flow  oscillation 
frequency  f  to  be  detected,  i.e.,  if  the  signal  is  almost  continuous.  One 
can  simply  treat  the  signal  as  a  continuous  hot-wire  anemometer  signal  to  obtain 
the  averages 


T 

U  =  y  f  U  (t)dt 
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T  J  (U( t)  -  U)n  dt 

o 


(1) 

(2) 


where  n  =  2,  3,  4  _  When  the  time  intervals  between  successive  signal 

bursts  are  long  compared  to  1 /fmax  (high  signal  dropout  rate)  and  are  unequal, 
these  equations  should  also  be  used  in  the  fashion  explained  below. 

First,  let  us  look  at  the  commonly  used  method  of  particle  averaging  for 
individual  particle  velocity  measurements.  The  averages  are  made  over  the 
number  of  signal  bursts  N  obtained  during  the  time  period  T: 

UN.  I=1  (3) 

N 

"7-  t,  (Bi-UN>n  (4) 

N 

where  n  =  2,  3,  4  _  These  averages  are  not  made  with  respect  to  time  and 

are  biased  unless  the  time  intervals  between  signal  bursts  are  equal.  McLaughlin 
and  Tiederman  (1973)  proposed  a  biasing  correction  that  is  based  upon  the  idea 
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that  higher  velocity  flow  carries  more  particles  through  the  focal  volume 
per  unit  time.  Thus,  more  high  velocity  signal  bursts  will  be  obtained  and 
will  be  too  high.  However,  high  velocity  particles  speed  less  time  in  the 
focal  volume  so  that  in  the  case  of  sampling  spectrum  analysis  signal  processing, 
the  chance  of  detecting  a  given  signal  burst  varies  as  ( +  Thus, 

this  effect  tends  to  cancel  the  above  mentioned  bias  for  particle  averaging. 

Durao  and  Whitelaw  (1975)  showed  that  if  the  Doppler  bursts  are  randomly  sampled 
before  particle  averaging,  the  bias  effects  are  reduced  significantly.  Even  so, 
particle  averaging  is  not  fundamentally  a  time  average. 

Consider  now  time-averaging  of  signals  according  to  equations  (1)  and 
(2),  even  though  the  signal  dropout  rate  may  be  large.  Only  ergodic  flows, 
whose  averaged  quantities  in  equations  (1)  and  (2)  become  independent  of  time 
for  large  T,  are  considered.  This  restriction  is  also  required  for  particle 
averaging.  The  last  sampled  signal  must  be  held  by  a  sample-and-hold  circuit 
until  a  new  signal  is  detected  for  time-averaging.  With  exception  of  the  instant 
at  which  a  new  signal  is  detected,  the  sampled-and-held  voltage  does  not  correspond 
to  the  actual  instantaneous  velocity.  However,  the  voltage  value  at  each  instant 
corresponds  to  the  instantaneous  velocity  at  some  instant  during  a  record  time  T 
for  an  ergodic  flow.  Since  any  averaging  process  removes  time  domain  dependency, 
it  does  not  matter  when  during  the  time  period  T  that  it  is  averaged.  It  is  unlikely 
that  a  given  signal  voltage  will  be  averaged  too  long  (Simpson  and  Chew,  1979). 

This  method  of  averaging  eliminates  the  need  for  the  high  velocity  flow  bias 
correction. 

The  mechanics  of  evaluating  a  true  time  average  in  this  research  made  use 
of  a  velocity  probability  histogram,  such  as  shown  in  Figure  5.  A  SAICOR  Model 
41  Correlator  and  Probability  Analyzer  was  used. 
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where  n  =  2,  3,  4  _  The  histogram  P(y)  is  constructed  by  sampling  the  U{ t) 

sample-and-held  signal  at  equal  intervals  in  time  At  for  the  period  T.  Thus 

the  histogram  reflects  a  true  time  integral  and  the  results  from  equations 

(5-7)  will  be  equivalent  to  those  from  equations  (1-2).  The  time  interval  At 

between  digital  samples  should  be  no  larger  than  the  shortest  time  between 

-4 

signal  bursts,  otherwise  some  data  will  be  lost.  For  example.  At  =  10  sec 
for  about  400  new  signals  per  second.  The  averaging  time  T  was  at  least  a  half 
minute,  so  at  least  12,000  new  data  signals  and  3x1 0  equal  time  interval 
samples  were  involved  for  one  histogram.  An  added  benefit  of  the  histogram 
approach  is  that  noise  can  be  detected  while  ?{u)  is  being  constructed.  If 
one  has  an  oscilloscope  display,  the  noise  will  cause  the  base  level  of  P (u) 
to  grow.  Thus,  the  resulting  P(U)  can  be  corrected  for  noise  or  the  discrim¬ 
inator  level  in  the  signal  processor  can  be  adjusted  on-line  and  a  new  P(y) 
constructed.  The  histograms  were  stored  on  digital  tape  and  analyzed  by  a  digital 
computer. 

These  results  are  not  believed  to  suffer  strongly  from  bias  errors.  First, 
there  is  no  bias  in  the  duration  of  a  detected  signal  due  to  the  flow  velocity. 

In  other  words,  the  time  that  the  highest  velocity  particle  spends  in  the  focal 
volume  is  always  large  enough  to  produce  a  sufficiently  large  vertical  voltage 
output  from  the  spectrum  analyzer.  Minimal  particle  seeding  was  used  for  the 
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Figure  5.  Typical  velocity  probability  histogram:  U  =  29.1fps,  u  =  27.0(fps) 

=  0.047,  =  2.66.  Discrete  velocity  bins  due  to  probability 

analyzer. 


Figure  6.  Freestream  velocity  distribution  along  the  tunnel  centerline. 
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best  SNR  and  data  sample  rate,  so  significant  finite  transit  time 
broadening  is  unlikely. 

Velocity  gradient  broadening  is  not  significant  for  any  data  presented 
here  (Simpson,  1976).  The  focal  volume  diameter  0.012  inches  and  length  0.140 
inches  are  small  compared  to  the  boundary  layer  thickness.  In  addition,  signals 
from  the  center  of  the  focal  volume  are  the  most  likely  since  the  scattered 
signals  are  the  most  intense.  Large-scaled  motions,  which  scale  on  the  boundary 
layer  thickness,  appear  to  dominate  the  structure  of  highly  turbulent  flows,  so 
strong  instantaneous  spatial  velocity  variations  are  unlikely.  In  any  event  as 
shown  below,  these  results  compare  favorably  with  hot-wire  anemometer  data 
obtained  in  regions  that  do  not  contain  significant  time  variation  of  the  flow 
direction. 

Since  no  spectra  has  previously  been  measured  in  the  separated  flow  zone, 
low  frequency  spectra  of  u  from  the  laser  anemometer  were  measured  using  a 
Princeton  Applied  Research  Model  4512  Fast-Fourier-Transform  Spectrum  Analyzer. 

II.  DESCRIPTION  OF  THE  TEST  FLOW 

All  data  were  obtained  at  atmospheric  pressure  and  77  +  %eF  flow  conditions. 
Figure  6  shows  the  free-stream  velocity  distributions  obtained  along  the  tunnel 
centerline  using  the  single-wire  probe.  This  distribution  was  repeatable  within 
2.9%  over  the  duration  of  these  experiments,  which  is  only  a  little  greater  than 
uncertainty  in  measuring  the  mean  velocity  with  a  hot-wire  anemometer  (+  2.4%). 
Figure  7  shows  the  non-dimensional  pressure  gradient  dCp/dx  along  the  centerline 
of  the  test  wall.  Here  Cp  =  2(p  -  p.)/pU^  =  1  -  (U^U^)2,  where  i  denotes  the 
free-stream  entrance  conditions  at  x  =  3".  A  five-point  local  least-squares  curve 

fit  of  data  was  used  at  each  streamwise  location  to  determine  this  derivative. 

P 
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Just  downstream  of  the  location  of  the  second  wall  jet  boundary  layer 
control  unit  (100  inches),  the  slope  of  the  static  pressure  gradient 
changes  sign.  Near  145  inches  the  pressure  gradient  drops  to  an  approxi¬ 
mately  constant  value  downstream. 

To  examine  the  two-dimensionality  of  the  mean  boundary  layer  flow, 
smoke  was  introduced  only  in  a  spanwise  portion  of  the  test  wall  boundary 
layer  at  a  given  time.  A  sheet  of  laser  light  produced  by  a  cylindrical 
lens  was  used  to  illuminate  the  smoke  across  the  tunnel.  Upstream  of  separa¬ 
tion,  negligible  spanwise  diffusion  of  the  smoke  was  observed,  indicating  no 
gross  flow  three-dimensionality.  Mean  velocity  profiles  at  several  spanwise 
locations  indicated  that  the  mean  velocity  was  two-dimensional  within  1%. 

Downstream  of  separation  greater  spanwise  diffusion  occurred,  so  that  downstream  of 
170  inches  no  nominally  two-dimensional  flow  remained.  On  the  basis  of  these 
observations,  the  wall  jet  and  suction  boundary  layer  controls  were  adjusted 
to  produce  a  nearly  two-dimensional  flow  pattern  downstream  of  separation. 

Smoke  flow  patterns  in  the  sidewall  and  comer  flows  were  symmetric  about 
the  channel  centerline. 

After  laser  and  hot-wire  anemometer  data  were  available,  examination  of 
the  two-dimensionality  was  done  by  evaluating  the  terms  in  the  two-dimensional 
continuity  equation  and  the  momentum  integral  equation.  In  the  first  method, 
the  differential  continuity  equation  was  written  as 

*  -  1  +  <UX  ♦  AX  -  V  /(I)'* 

R  was  computed  by  finding  the  gradient  of  V  with  respect  to  y  and  also  the 
change  in  U  with  respect  to  X  at  a  constant  y  location.  Only  where  LDV  data 
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were  available  was  this  method  useful,  because  it  requires  good  V  data. 
For  many  streamwise  locations  R  lies  between  +0.5  and  -0.5  in  the 


middle  region  of  the  boundary  layer,  with  an  uncertainty  of  +  0.46. 

Nearer  the  wall,  AU  becomes  relatively  more  uncertain  while  3V/3y  is 
more  uncertain  in  the  outer  region.  As  a  result  one  can  expect  greater 
uncertainties  in  R  in  these  regions.  Thus,  at  least  in  the  middle  region  of 
the  boundary  layer,  the  flow  is  two-dimensional  within  the  uncertainty  of 
evaluating  R.  This  is  a  stringent  test  since  it  is  based  only  on  the  local 
flow  field,  but  it  suffers  from  the  disadvantage  of  needing  to  differentiate 
experimental  velocity  distributions. 

On  the  other  hand,  the  momentum  integral  equation  provides  a  global 
test  based  on  conservation  of  momentum  over  a  large  flow  volume.  In  this 
method,  the  momentum  integral  equation  was  again  integrated  in  the  x-di recti  on 
to  yield 


dlL 


1  -a;  dx  -  / 


x 

d  x  +  / 


x0 


J  _3  (u ' 2  -v'2) 
0  3x 


dx 


(8) 


where  the  l.h.s.  is  PL  and  the  r.h.s.  is  PR.  The  last  term  of  PR  is  due  to  the 
normal  stresses  and  its  effect  in  the  vicinity  of  separation  has  been  shown 
(Simpson  et  al.,  1977)  to  be  significant.  Using  the  experimental  data,  PL  and 
PR  of  eqn.  (8)  were  evaluated  with  and  without  the  normal  stress  term  and  the 
ratio,  j^-  -  1,  was  computed  for  both  the  cases.  The  distributions  on  the  ratio 
in  the  streamwise  direction  are  shown  in  Fig.  8  along  with  estimated  uncer¬ 
tainty  bands.  The  ratio  computed  without  the  normal  stress  term  is  within 
+  0.16  up  to  122  inches,  indicating  that  the  flow  is  reasonably  two-dimensional. 
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Figure  8.  Fractional  imbalance  of  the  momentum  integral  equation, 
PL/PR  -  1,  vs.  X.  Associated  uncertainty  denoted  by 
vertical  band. 
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Downstream  the  normal  stresses  play  an  more  important  role,  although 
they  are  not  large  enough  to  account  for  the  imbalance  far  downstream 
of  separation. 

IV.  EXPERIMENTAL  RESULTS  FOR  THE  MEAN  FLOW 
IV. 1  Mean  Velocity  Profiles 

Mean  velocity  profiles  were  obtained  in  the  unseparated  upstream 
boundary  layer  and  the  outer  part  of  the  separated  flow  using  single  wire 
and  cross-wire  hot-wire  anemometer  probes.  The  di recti onally-sensitive 
laser  anemometer  provided  velocity  profiles  in  the  separated  zone  and  the 
region  immediately  upstream. 

Figures  9a  and  9b  show  the  streamwise  mean  velocity  profiles  for 
steady  flow  for  a  few  typical  stations  in  the  near  separation  and  the  sepa¬ 
rated  regions  obtained  using  all  three  different  techniques.  There  is  good 
agreement  among  these  measurements,  with  the  maximum  discrepancy  among  them 
being  about  6  to  7%.  In  the  separated  region  only  the  laser  anemometer 
measurements  are  meaningful.  Table  1  presents  the  experimental  uncertainties 
for  each  measured  quantity  as  determined  by  the  method  of  Kline  and  McClintock 
(1953).  As  shown  by  Simpson  and  Chew  (1979),  the  laser  anemometer  results 
obtained  on  different  days  at  the  same  location  indicate  a  high  level  of  data 
repeatability. 

Figures  10  and  11  show  non-dimensional  streamwise  mean  velocity  profiles 
across  the  boundary  layer  at  various  streamwise  stations  in  linear  and  semi- 
logarithmic  co-ordinates.  These  results  were  obtained  by  smoothing  a  curve 
between  the  laser  and  valid  cross-wire  data.  While  the  smoothing  was  a  some¬ 
what  subjective  procedure,  one  can  see  from  Figures  9a  and  9b  that  this 
procedure  basically  just  eliminated  a  few  scattered  data  points. 
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Figure  9  (a).  Non-dimensional  streamwise  mean  velocity  profile  data 
from  the  laser  and  hot-wire  anemometers.  Note  dis¬ 
placed  ordinates. 
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Figure  9  (b).  Non-dimensional  streamwise  mean  velocity  profile 
data  from  the  laser,  single  hot-wire,  and  cross 
hot-wire  anemometers.  Note  displaced  ordinates. 
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Smoothed  results  from  laser  and  hot-wire 
anemometer  data.  Note  displaced  ordinates 
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Figure  11(a).  Non-dimensional  streamwise  mean  velocity 
"profiles:  single  hot-wire  results. 
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Figure  11(b).  Non-dimensional  streamwise  mean  velocity 
profiles:  single  hot-wire  results. 
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11(c).  Non-dimensional  streamwise  mean  velocity 
profiles:  smoothed  laser  and  valid  single 
hot-wire  anemometer  results. 


XLIncmx] 


X 

x 


Y/IO.TR 


Figure  11(d).  Non-dimensional  streamwise  mean  velocity 
profiles:  smoothed  laser  and  valid  single 
hot-wire  anemometer  results. 


Figure  12  shows  non-dimensional  laser  anemometer  and  cross-wire 
anemometer  results  for  the  normal  velocity  component  V  just  upstream  of 
separation  and  in  the  separated  region.  At  most  streamwise  locations, 
there  is  good  agreement.  However,  as  shown  in  Table  1,  there  is  a  fairly 
large  uncertainty  in  the  cross-wire  result,  mainly  because  of  the  uncer¬ 
tainty  of  the  probe  orientation  with  respect  to  the  test  wall.  Therefore, 
the  laser  anemometer  results  are  more  reliable. 

IV. 2  Turbulence  Quantities 

_  2 

Figures  13,  14,  and  15  show  u'/U^,  v'/U^,  and  -uv/U^  vs.  y/6,  respec¬ 
tively.  The  agreement  between  the  laser  and  cross-wire  anemometer  results  is 

good  with  the  apparent  discrepancies  being  due  to  the  experimental  uncertain- 

_  2 

ties  shown  in  Table  1.  The  discrepancies  in  the  -uv/U^  plots  are  the  great¬ 
est  due  to  the  uncertainty  in  orientation  of  the  cross-wire  probe  with 

~2  _ 

respect  to  the  test  wall.  Since  u  and  v  are  much  larger  than  -uv,  only 

a  very  small  misalignment  is  required  to  produce  a  much  different  -uv  result. 

Figures  16  (a-e)  show  profiles  obtained  by  smoothing  a  curve  between  the  laser 

anemometer  and  valid  cross-wire  data.  Figure  17  shows  u'/U^,  v'/U^,  and 

-uv/U^  profiles  upstream  of  the  near  separation  zone.  Figures  18  (a)  and  (b) 

show  u'/U^  obtained  from  the  single  wire  anemometer. 

IV. 3  Upstream-downstream  Intermittency 

Only  the  directionally-sensitive  laser  anemometer  results  from  these 

measurements  give  information  on  the  fraction  of  time  that  the  flow  moves 

downstream  or  y  .  This  quantity  is  the  fraction  of  the  area  of  the  velocity 
pu 

probability  histogram  that  has  a  positive  velocity.  The  directionally-insen- 
sitive  hot-wire  anemometer  cannot  yield  ypu  values  (Simpson,  1976). 
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Normal  hot-wire: 


Cross  hot-wire: 
(including  misalignment 
uncertainty) 

Laser  anemometer: 


Position  from  wall: 

Skin  friction 
coefficient  C^.: 


U  -  2.4%,  7-7% 

U  -  3.2%,  u2  -  10% 

v2  -  11%,  -U7  -  17%,  5y  -  0.1 
U,  V  -  0.2  fps;  7  and  v2 

+  4%  max.  profile  value;  uv  -  6% 
max.  profile  value;  skewness  -  0.1; 
flatness  -  0.2;  YpU  -  0.1  exp(-U2/27) 
Ypv  -  0.1  exp(-V2/2v2) 

-  0.002  inches 

Ludwieg-Tillman,  -  6.5%; 

Preston  tube,  -  8.5%; 

Surface  hot-wire,  -  12%. 


Table  1.  Estimated  uncertainties  of  measured  quantities. 
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Figure  12(b).  Mean  velocity  V/U^  vs.  y/6  profiles,  laser 
and  cross  hot-wire  anemometer  data. 
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Figure  13(a).  Streamwise  turbulence  intensity  u'/U^  vs.  y/6 

profiles:  laser,  single  hot-wire,  and  cross  hot¬ 
wire  anemometer  results.  Note  displaced  ordinates. 
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Figure  13(b).  Streamwise  turbulence  intensity  u'/U  vs.  y/6 
profiles:  laser,  single  hot-wire,  and  cross 
hot-wire  anemometer  results.  Note  displaced 
ordinates . 
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Figure  14(b).  v'/U^,  vs.  y/6  profiles:  laser  and  cross 

hot-wire  anemometer  results.  Note  displaced 
results. 
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Figure  15(a). 


_  2 

Reynolds  shearing  stress  -uv/U^  vs. 
profiles:  laser  and  cross  hot-wire 
anemometer  results.  Note  displaced 


y/6 

ordinates. 
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Figure  15{b).  Reynolds  shearing  stress  -uv/U^  vs.  y/6 

profiles:  laser  and  cross  hot-wire  anemometer 
results.  Note  displaced  ordinates. 
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Figure  16(b).  Smoothed  u'/U  vs.  y/6  results  from  hot-wire 

and  laser  anemometers.  Note  displaced  ordinates. 
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Figure  16(d).  Smoothed  v'/U^  vs.  y/6  results  from  hot-wire 

and  laser  anemometers.  Note  displaced  ordinates 
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Figure  16(e),  Smoothed  Reynolds  shearing  stress  -uv/U^ 
vs.  y/5  results  from  laser  and  hot-wire” 
anemometers.  Note  displaced  ordinates. 
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Figure  16(f).  Smoothed  Reynolds  shearing  stress  -uv/U^ 

vs.  y/5  results  from  laser  and  hot-wi  re°°ane- 
mometers.  Note  displaced  ordinates. 
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Figure  17(a). 


Cross  hot-wire  anemometer  results  well  upstream 
of  separation:  u’/U^. 
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Figure  17(b).  Cross  hot-wire  anemometer  results  well 
upstream  of  separation:  v’/U^. 
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Figure  18(b).  u'/U  vs.  y/5  single  hot-wire  anemometer 
results. 
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Figure  19  shows  the  distributions  of  the  intermittency  across  the 
boundary  layer  for  the  region  approaching  separation  and  downstream  of  it. 

The  reverse  flow  first  starts  appearing  at  122.6  inches  but  becomes  clearly 
observable  beyond  127  inches.  Further  downstream,  the  backflow  intensifies 
and  also  spreads  outwards  from  the  wall,  y  reaches  the  lowest  value  of 
approximately  0.05  at  the  last  station  of  measurement  in  the  separated 
region,  where  backflow  extends  up  to  about  60%  of  the  boundary  layer  thick¬ 
ness.  The  distributions  in  the  separated  region  are  trough-shaped  near  the 
wall  showing  that  the  maximum  amount  of  reverse  flow  occurs  slightly  away 
from  the  wall.  This  is  consistent  with  the  velocity  profile  shape  that  shows 
that  the  highest  velocity  for  the  backflow  is  reached  at  a  point  slightly 
away  from  the  wall.  However,  as  shown  in  Table  1,  the  uncertainty  in  ypu 
becomes  large  as  the  mean  velocity  approaches  zero,  so  one  cannot  place  too 
much  emphasis  on  this  coincidence.  Figure  20  shows  the  decay  of  ypu  near 

the  wall,  y  ,  as  a  function  of  the  streamwise  co-ordinate.  As  mentioned 
puo 

earlier  the  reverse  flow  is  first  measured  at  122.6  inches  and  thereafter 
its  persistence  increases  continuously,  until  it  reaches  a  level  where  it 
exists  90%  of  the  time  after  which  its  rate  of  increase  diminishes. 

IV.  4  Higher-order  Turbulence  Correlations 

To  investigate  the  effect  of  separation  on  higher  order  structure 
functions,  the  third  and  fourth  moments  given  by  equation  (7)  were  calcu¬ 
lated  from  the  u  and  v  LDV  histograms.  Simpson  and  Chew  (1979)  showed  that 

—k  —k  3/2  -3  3/2 

the  skewness  factors,  Su  =  (uJ)/(u  )  and  Sy  =  (vJ)/(v  )  ,  and  flatness 

factors,  Fu  =  (u*)/(?)2  and  Fy  =  (u^)/(v^  were  about  +  0.1  and  +  0.2 
uncertain.  Data  obtained  on  different  days  were  in  close  agreement,  with  the 


54 


13US75 

13%iUD 

13E.7SD 

1**.STS 

1S637S 

17HBT5 


19(b). 


Streamwise  fraction  of  time  flow  moves 
downstream,  ynu  vs.  y/6. 

ru 


56 


scatter  being  within  these  uncertainty  levels. 

~2~ 

For  purposes  of  comparison  and  for  additional  information,  u  v  and 
v  triple  correlation  data  were  obtained  from  the  cross-wire  anemometer. 

The  main  source  of  uncertainty  in  these  measurements  is  the  drift  of  the 
mean  voltage  level  in  the  multipliers.  This  was  kept  to  a  minimum  by 
adjusting  the  offset  voltage  before,  several  times  during,  and  after  taking 
a  set  of  data,  so  that  a  zero  voltage  input  produced  a  zero  voltage  output. 
During  data  reduction  a  correction  was  applied  for  the  offset  voltage. 

Table  2  gives  the  maximum  uncertainties  for  each  streamwise  location  that 
data  are  presented.  All  data  which  were  greater  25%  uncertain  are  not 
presented  here.  This  arbitrary  uncertainty  limit  is  not  really  very  high, 
considering  that  third-order  correlations  are  expected  to  have  high  uncertain¬ 
ties. 

Figure  21  shows  the  skewness  factor  Sy  results  obtained  by  laser  and 
cross-wire  anemometers  for  several  streamwise  locations.  The  agreement  between 
the  two  types  of  experimental  results  away  from  the  wall  at  a  given  location  is 
generally  within  the  estimated  uncertainties.  In  the  separated  zone  the  hot¬ 
wire  measurements  were  confined  to  the  outer  region  where  the  instantaneous  flow 
direction  differed  less  than  45°  from  the  mean  flow  direction. 

IV. 5  Skin-friction  Results 

Three  different  ways  of  deducing  the  near  wall  shearly-stress  distri¬ 
bution  were  used:  the  Ludwi eg -Tillman  skin-friction  correlation,  a  Preston 
tube,  and  the  surface  hot-wire  gage  described  in  section  II. 4  above.  Figure 
22  shows  the  results  from  these  three  methods,  which  are  in  agreement  within 
the  uncertainties  given  in  Table  1.  Table  3  gives  the  Ludwi eg-Ti 11 man  results. 
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Table  2:  Uncertainties  for  the  diffusion  results,  y/6  positions  given  in 
parentheses . 


X 

Inches 

value  at  the  loca¬ 
tion  where  the  un¬ 
certainty  was 
computed 

Maximum  absolute 
value 

Maximum  value  of 
the  absolute 
uncertainty 

u2v  , 
(ft/sec) J 

3 

v  , 

(ft/sec) J 

u2v  3 
(ft/sec) J 

7 

(ft/sec)3 

u2v  ~ 

( ft/sec r 

v3  , 
(ft/sec)' 

_ 

(ftTsec) 

31.688 

0.58 

(0.401) 

3.93 

(0.465) 

4.78 

6.43 

0.2 

0.15 

62.36 

86.5 

1.07 

(1.221) 

1.03 

(1.396) 

14.53 

10.09 

0.91 

1.0 

66.73 

117.625 

3.09 

(0.27) 

3.02 

(0.407) 

11.49 

3.5 

1.7 

0.25 

51.86 

126.75 

12.72 

(0.623) 

4.83 

(0.837) 

16.6 

8.59 

1.02 

0.33 

48.94 

131.0 

13.38 

(0.714) 

3.44 

(0.88) 

13.86 

4.69 

1.37 

0.5 

47.97 

144.0 

9.67 

(0.895) 

4.13 

(0.933) 

11.24 

6.2 

1.03 

0.39 

47.06 

156.375 

6.97 

(0.966) 

1.53 

(1.073) 

14.27 

9.74 

1.57 

0.88 

45.33 
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Friction  factor  Cf/2  vs,  x  from  the  several  experimental 
methods:  ■,  Ludwig  -Tillman  lawjd,  surface  hot-wire  gage; 
♦  ,  Preston  tube. 


CM 

1 

1  o 

o 

o 

o 

o 

o 

o 

o 

o 

<\i  a 

O' 

o 

OC 

o 

4- 

4! 

in 

o 

«-• 

© 

«4» 

•  in 

r-» 

^4 

m 

CM 

< 

cn 

o 

O' 

CC  sD 

cn 

f-4 

00 

a> 

4 

4 

CM 

o 

o 

C' 

1  m 

o 

cc 

in 

•£> 

CM 

o 

m 

in 

o-  >»■ 

CM 

o 

4* 

oc 

4> 

rH 

O 

o 

o 

X 

1  o 

o 

o 

O' 

a 

r- 

m 

m 

or 

nJ- 

4- 

cn 

rH 

o 

o 

o 

o 

o 

o 

"o 

•  <M 

CM 

CM 

H 

rH 

r4 

^4 

r4 

i 

i 

1 

1  CL 

r» 

© 

cn 

cn 

N 

cn 

CO 

tn 

cn  O' 

o 

CM 

vC 

m 

cn 

© 

QD 

4 

© 

m 

g 

1  r- 

o 

m 

CM 

*1* 

m 

r4 

^4 

CM 

cn  cc 

rH 

r*4 

cr 

r 

r- 

cn 

rH 

o 

QD 

m 

1  cn 

fr> 

co 

cn 

cn 

cn 

4 

<4* 

>o 

f"  O'  o 

CM 

o 

in  o 

rH 

H 

4 

<0 

4 

X 

1  • 

• 

• 

• 

• 

• 

• 

• 

• 

•  • 

• 

• 

m 

0 

• 

• 

• 

• 

• 

• 

I  r4 

1 

1 

rH 

rH 

rH 

H 

H 

*■4 

»—4 

*-4  i-H 

CM 

CM 

CM 

CM 

CM 

cn 

4 

in 

© 

1 

1 

1  <\J 

CM 

o 

a 

O ' 

4> 

cn 

O' 

cn  n- 

00 

r- 

O' 

© 

CM 

r4 

O' 

r- 

4 

CD 

u>* 

I  f— < 

o 

n- 

m 

o 

o 

O' 

rH 

sO 

cd  >r 

4- 

4- 

h- 

00 

00 

<o 

in 

4 

4 

a: 

1  CD 

r- 

4 

cr 

o 

•c 

o 

r4 

oc  m 

CM 

OD 

4* 

oc 

cn 

cn 

OD 

© 

CM 

•o 

1  cn 

m 

<c 

cn 

4 

© 

<c 

CM 

4) 

co  r* 

O' 

<c 

o 

rH 

rH 

CM 

o 

4) 

i 

rH 

rH 

«v 

CM 

CM 

cn 

in 

00 

O'  ^ 

CM 

cn 

4* 

© 

<£) 

© 

h* 

■a 

0D 

00 

i 

i 

^4 

rH 

^4 

rH 

rH 

^4 

^4 

rH 

rH 

rH 

H 

• 

1  © 

QD 

r- 

CD 

rH 

fO 

rH 

CM 

nj 

M  (M 

CM 

rH 

4- 

cn 

© 

r* 

© 

rH 

o 

• 

1  © 

CM 

o 

O' 

o 

o 

0D 

00 

H 

O  M 

O' 

o 

4* 

in 

4 

O' 

© 

o 

in 

CM 

to 

1  >D 

f>- 

CD 

© 

rJ 

4 

-<  «c 

cn  ** 

m 

rH 

4* 

sD 

41 

4 

o 

in 

rH 

<-o'E 

1  o 
|  • 

O 

• 

o 

• 

o 

• 

o 

• 

H 

• 

rH 

• 

CM 

• 

4* 

• 

>c  a 

•  • 

o 

• 

CM 

• 

>0 

• 

o 

• 

O 

• 

4) 

• 

rH 

• 

o 

• 

CM 

• 

cn 

• 

w’ 

i 

i 

rH 

iH 

CM 

CM 

CM 

cn 

4 

■c 

CD 

i 

• 

u  • 

1  r- 

in 

cn 

rH 

CM 

O' 

o 

^4 

CD  CM 

cn 

CM 

4* 

CM 

vD 

in 

4 

0D 

■a 

H 

1  f-t 

H 

<£) 

o 

<o 

in 

O' 

O' 

od  >r 

CM 

4* 

CM 

n- 

© 

cn 

■o 

~o 

4 

,  *-z 

1  in 

nO 

in 

© 

r>- 

»  O' 

CM  CO 

cn  oc  O* 

<o 

0D 

© 

o 

OD  ® 

t ^ 

o 

00 

V-0  M 

1  • 

• 

• 

• 

• 

• 

• 

• 

• 

•  • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

w 

i 

H 

*-4 

2 

2 

2 

cn 

cn 

4- 

© 

m 

in 

>o 

o 

CM 

i 

i 

rH 

rH 

i 

i 

0 

1  CM 

4 

H 

4> 

CM 

CM 

O' 

m 

CD 

CC  O' 

•o 

o 

4" 

OD 

o 

©  © 

CO 

%o 

0D 

i  © 

h- 

m  n- 

O' 

o 

cn 

CM 

m  O' 

4* 

oc 

o 

rH 

CM 

r- 

o 

4 

rH 

.  ©z 

1  4 

in 

in 

O' 

CM 

OD 

cn  h- 

o 

CM 

OD 

OD 

<o 

QD 

r- 

r- 

m 

in 

*-©  *h 

|  • 

• 

• 

0 

• 

• 

0 

• 

• 

•  • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

i 

rH 

T-* 

CM  CM 

CM 

cn 

cn 

4* 

4 

m 

m 

4) 

© 

CM 

i 

i 

r4 

X  — * 

i 

i 

o  • 

i 

V  U 

1  cd 

CD 

o 

CD 

o 

CM 

vC 

0* 

CD  CM 

in 

H 

O' 

■£> 

in 

in 

4 

OD 

1-4 

LL  UJ 

I  m 

4 

m 

m 

4 

m 

m 

in 

cn 

r- 

4- 

rH 

4- 

ao 

in 

o 

<o 

CM 

CM 

Z  l/> 

1  • 

• 

• 

• 

• 

• 

• 

0 

• 

•  • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

»-< 

1  in 

o 

o 

4 

CM 

H 

CVJ 

cn 

in  4- 

4- 

4* 

cn 

cn 

CM 

CM 

cv 

rH 

•H 

rH 

13  -I 

i 

i 

1 

• 

* 

i  i 

i 

| 

i 

i 

i 

i 

i 

i 

i 

i 

Q  — 

• 

i 

• 

i 

« 

• 

o 

1  o 

o 

o 

o 

o 

o 

o 

o 

o 

o  o 

o 

o 

© 

o 

o 

o 

o 

3 

o 

o 

u.  UJ 

1  o 

o 

o 

o 

o 

o 

o 

o 

o 

o  Oin 

o 

o 

in 

o 

o 

in 

O 

o 

o 

Z  to 

1  •£> 

00 

CM 

n 

<£ 

o 

CD 

cn 

O  -4 

H 

cn 

in 

r4 

4 

4 

m 

-o 

r*- 

•■4  V* 

|  • 

• 

• 

• 

• 

• 

• 

• 

• 

•  • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

x  . 

1  CM 

r- 

m 

or- 

fH 

o 

<c 

a 

ir  cn 

CM 

^H 

o 

CC 

OC 

4) 

in 

4 

cn 

1  IT. 

IT 

4. 

r- 

N 

in 

in  ir 

in 

in 

4* 

4 

4 

4 

4 

4 

4 

4 

a 

i 

i 

i 

1 

1  m 

o 

c 

o 

o 

C 

©  o 

o 

l Tv  O 

c 

in 

in 

m 

in 

©  o 

in 

in 

in 

• 

•  N 

in 

tr 

o 

o 

IT 

in  o 

^-4 

h-  C 

c 

CM 

CM 

N 

r- 

c 

in 

r 

to 

1  r 

cv 

o 

IT 

CV 

CM 

IT 

cn 

m  m 

m 

•c 

r*~ 

fr 

CC 

IT 

r- 

tt 

CD 

X  Z 

|  0 

• 

• 

• 

0 

• 

• 

0 

• 

•  • 

• 

• 

• 

0 

• 

• 

• 

• 

• 

• 

1  H 

rH 

rH 

4 

CM 

4 

4 

r- 

CM  CD 

o 

<\j  N 

o 

r4 

4 

CC 

4 

4> 

c 

w 

1  r4 

CV 

cn 

4 

IT 

<c 

r- 

er 

o 

I-4  i-4 

CM 

CM 

CM 

CM 

cn 

cn 

cn 

4 

in 

i 

H 

^4  H 

r4 

rH 

^4 

r4 

r4 

H 

rH 

rH 

^H 

62 


1 


Table  3:  Parameters  of  the  mean  flow  development 
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The  Preston  tube  and  Ludwieg-Tillman  methods  require  the  existence 
of  a  universal  logarithmic  law-of-the-wall  velocity  profile.  The  data 
obtained  using  the  surface  hot-wire  gage  are  not  dependent  on  the  require¬ 
ment  of  a  logarithmic  wall  region.  This  suggests  that  the  law  of  the  wall 
is  valid  until  the  location  where  YpU  is  first  less  than  one  near  the  wall. 
These  results  are  in  agreement  with  results  of  Simpson  et  al . ,  (1977). 

IV. 6.  Data  Tabulation 

These  data  are  tabulated  in  Table  3  and  in  the  Appendix.  These  data 
are  recorded  on  magnetic  tape  in  the  format  required  for  the  1980-81  AF0SR- 
HTTM-Stanford  Conferences  on  Complex  Turbulent  Flows,  a  copy  of  which  is  on 
file  in  the  Thermosciences  Division  of  the  Stanford  University  Department  of 
Mechanical  Engineering. 


V.  DISCUSSION 

V.l  Mean  Velocity  Distribution 

Figure  23  shows  the  mean  streamline  pattern  for  the  flow  in  the  vicinity 
of  separation.  Note  that  in  the  backflow  region  the  turbulence  level  is  very 
high  compared  to  the  mean  flow,  so  these  mean  streamlines  do  not  represent  the 
average  pathlines  for  elements  of  fluid.  As  discussed  in  section  VI  below,  it 
appears  that  the  fluid  in  the  backflow  does  not  come  from  far  downstream  as  the 
streamlines  may  suggest,  but  is  supplied  fairly  locally. 

Figures  24  (a)  and  (b)  show  that  the  U+  vs  y+  law-of-the-wall  velocity 
profile  holds  all  along  the  flow  channel  when  the  Ludwieg-Tillman  skin  friction 
values  are  used.  Although  no  wall  proximity  corrections  to  the  hot-wire  data 
were  applied  in  the  viscous  sublayer,  the  U+  s  y+  relationship  is  obeyed  rather 
well.  Oka  and  Kostic  (1972)  noted  that  hot-wire  measurements  are  only  influ¬ 
enced  by  flow  interference  and  conduction  to  the  test  wall  for  y+  <  4,  which 
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explains  why  the  present  data  for  y+  >  4  obey  the  viscous  sublayer 
equation  so  well.  Upstream  of  the  vicinity  of  separation,  the  usual 
logarithmic  form  for  y+  >  30  holds 

U  =  In  y  +  5.0  (9) 

as  shown  in  Figures  24. 

Perry  and  Schofield  (1973)  proposed  universal  empirical  correla¬ 
tions  for  the  inner  and  outer  regions  of  adverse  pressure  gradient  boundary 
layers  near  separation.  Their  correlations  apply  to  all  types  of  adverse 
pressure  gradient  boundary  layers  irrespective  of  whether  they  are  in 
equilibrium  or  not,  but  with  the  restriction  that  the  ratio  ("uv)max/,|JT 
must  exceed  1.5.  They  proposed  the  defect  law  for  the  outer  flow  as 


IL  -  U 


f 2(n?).  where  *  y/A 


Uo» 

.  _  00  I 
A  "  us  C 


C  is  a  universal  constant  given  by  C  =  /  fg^g)  and  found  empirically 

to  be  equal  to  2.86.  The  inner  law  was  defined  as 

2 

IJ  ‘  h  *  W*  nl  =  e  ’  e  = 

T  UMP 

2 

where  h  is  a  constant  and  U^p  and  L  are  described  later. 

The  condition  for  the  overlap  between  the  inner  and  the  outer  region 

lead  to  the  following  relations 

/11  \3/2  /  \l/2  /•  1  \ 


/u  Y/d  (  \/d  /u  \ 
U_e  0  47  (-1)  (jL]  J  -1} 

0A7VJ  V6!/  \U»/ 
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and 


1/2 

f-j  (r^)  =  6.4 


(14) 


(15) 


Equation  (13)  was  used  to  obtain  the  values  of  the  velocity  scale  factor  U 


by  drawing  a  Clauser-type  chart  for 


k-(kT 


with  U$  as  the  parameter. 


All  the  parameters  obtained  for  Perry's  correlations  are  given  in  Table  4. 

_  2 

The  condition  (-uv)max/uT  >1*5  was  satisfied  by  the  data  for  the 
region  downstream  of  x  =  105  inches.  Hence,  Perry  and  Schofield's  correla¬ 
tions  were  tried  for  the  locations  downstream  of  105  inches  where  the  profiles 
of  mean  velocity  as  well  as  those  of  normal  and  shear  stresses  were  available. 
The  data  for  the  normal  stresses  are  also  required  since  Perry  and  Schofield 
neglected  the  normal  stresses  term  in  the  momentum  equation  while  computing 
the  shear  stress  profiles  from  the  mean  velocity  profiles.  It  was  shown  later 
by  Simpson  (1975)  that  the  normal  stresses  term  plays  a  significant  part  in 
both  the  momentum  and  the  turbulence  energy  equations  for  flows  approaching 
separation.  The  normal  stresses  effects  have  been  considered  in  a  way  as 
discussed  by  Simpson  et  al.  (1977)  and  in  accordance  the  pseudo-shear  stress 


UMp  is  defined  as 


2  00  "T  — x- 

UMp  =  -uv  +  /  9 Cu2  -  v  )  dy 

y  3x 


(16) 


L  J  max 

and  L  is  defined  as  the  distance  from  the  wall  to  the  maximum  in  the  pseudo¬ 
shear  stress  profile. 

Figures  25  and  26  show  the  velocity  profiles  plotted  in  the  inner  and 
outer  layer  co-ordinates.  The  inner  law  correlation  given  by  eqns.  (12)  and  (14) 
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Table  4.  Experimental  values  for  the  parameters  used  in  the  Perry  and 
Schofield  correlation. 


ill .125 


118.5 


126.75 


u/u 

s  00 

BiMiiiiii 

yc/A 

EZHI 

'  0.67 

15.86 

15.89 

0.0889 

1.077 

0.78 

16.80 

16.64 

0.0421 

1.103 

0.86 

21.33 

18.21 

0.0392 

1.183 

1.02 

22.47 

18.66 

0.0151 

1.212 

'  .99 


0.340  | 


0.284 


0.423 


0.617 
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slope  «  8.05 


Figure  25.  Perry  and  Schofield  inn^r  region  correlation  for  the  present 

data  near  separation,  U-  h  vs.  (y/e)’/S  equations  (12)  and  (14) 
given  by  solid  lines  for  6.4  and  8.05  slopes. 


y/a 


Figure  26.  Perry  and  Schofield  outer  region  correlation  for 
the  present  data  near  separation,  (U^  -  U)/U  vs. 
y/A,  eqn.(10);  solid  line  is  mean  line  from  Perry 
and  Schofield  (1973). 


seem  to  be  satisfied  reasonably  well,  although  the  higher  slope  of  8.05 
satisfies  the  upstream  most  profiles  better  (Simpson  et  al . ,  1977).  In 
the  region  near  the  wall,  eqn.  (12)  takes  the  usual  logarithmic  form  of 
eqn.  (9).  By  matching  the  logarithmic  and  the  half  power  regions.  Perry 
and  Schofield  obtained  the  expression  for  the  point  of  tangency  as  yc  =  0.58e. 

As  shown  in  Table  4  the  predicted  point  of  tangency  moves  toward  the  wall  as 
one  proceeds  downstream,  indicating  that  the  extent  of  the  logarithmic  region 
gradually  decreases,  which  can  also  be  seen  in  Figure  11c.  The  present  data 
satisfy  the  other  matching  condition  given  by  eqn.  (15)  to  a  reasonable  extent. 
The  present  data  upstream  of  intermittent  backflow  lie  within  the  band  repre¬ 
sented  by  the  scatter  in  the  data  plotted  by  Perry  and  Schofield. 

Following  Strickland  and  Simpson  (1973),  the  velocity  profiles  in  the 
separated  region  were  normalized  to  see  whether  the  profiles  in  the  outer  region 
show  any  resemblance  to  those  observed  in  mixing  layers.  For  this  purpose 
U_U0  was  plotted  as  a  function  of  as  shown  in  Figure  27  for  a  few 

V“o  0 

stations  downstream  of  127  inches.  In  the  case  of  the  mixing  layer,  y*  repre¬ 
sents  the  distance  from  the  center  of  the  mixing  layer  and  here  the  location 
where  the  Reynolds  Shear  Stress  -uv  reaches  a  maximum  was  considered  as  the 
center  of  the  shear  layer.  Xq  is  the  streamwise  distance  from  a  reference 
point  and  in  the  present  studies  x  =  88"  was  taken  as  the  reference  point.  Uq 
is  equal  to  twice  the  velocity  at  the  center  of  the  shear  layer  minus  the  free- 
stream  velocity  and  o'  is  a  constant.  Also  shown  is  the  curve  obtained  by 
Halleen  (1964)  for  a  mixing  layer.  An  error  function  type  of  distribution  repre¬ 
sented  by 

U  -  U„ 


U, 


1  +  erf 


(17) 
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Figure  27.  Present  data  in  mixing- layer- type  co-ordinates 
solid  line,  eqn.(17)  with  a'=  24  and  X  =  88"; 
dashed  line,  Halleen's  (1964)  mixing  layer  cor 
relation. 
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is  also  plotted  for  comparison.  There  is  good  agreement  of  the  data  with 
these  distributions,  confirming  the  analogy  with  mixing  layers.  However,  a1 
is  about  24  while  Halleen  obtained  a  value  of  about  17.5. 

As  one  can  see  in  Figure  lid,  there  is  some  profile  shape  similarity 
for  the  backflow  mean  velocity  downstream  of  138  inches.  Figure  28  shows  a 
good  correlation  when  normalized  on  the  maximum  negative  mean  velocity  and 
its  distance  from  the  wall  N.  A  slightly  poorer  correlation  results  when  6  is 
used  instead  of  N.  The  U+  vs.  y+  law-of-the-wall  velocity  profile  is  not  con¬ 
sistent  with  this  correlation  since  both  and  N  increase  with  streamwise 
distance,  while  the  law-of-the-wall  length  scale  v/UT  varies  inversely  with  its 
velocity  scale  UT.  The  data  of  Simpson  et  al .  (1977)  for  the  one  available 
location  are  also  shown  to  be  in  fair  agreement  with  this  correlation. 

An  attempt  was  made  to  see  if  the  mean  velocity  profiles  downstream  of 
separation  could  be  composed  of  the  "law-of-the-wake"  (Coles  and  Hirst,  1969) 

w(y/6)  =  2  sin^2fiy  C18) 


and  a  similarity  distribution  for  the  remaining  wall  flow.  Figure  29,  which  is 

o 

a  plot  of  U/U^  -  sin  (Try/26)  vs.  y/6,  shows  the  remainder  for  the  wall  flow. 
There  is  no  significant  profile  similarity. 

Another  attempt  was  made  to  scale  the  wake  function  by  using  the  maximum 
backflow  velocity  and  the  free-stream  velocity  before  subtracting  it  from  the 
velocity  profile.  This  was  done  as  follows: 


u  _ 

r  u~  +  un[ 

u 

u 

00 

00 

1 

I 


uj(y/6)  - 


+  R(y/«) 


(19) 


where  R(y/6)  can  be  called  a  "backflow"  function.  Furthermore,  another  function 
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B(y/<$)  =  R(y/6)  was  formed  so  that  B(y/6)  has  definite  limits  of  1  and 
kNi 

0  at  y/6  =  0  and.l,  respectively.  The  plots  of  these  function  R(y/6)  and 


B(y/6)  are  shown  in  Figures  30  and  31.  They  neither  show  any  similarity  nor 
small  values  in  the  outer  region.  This  leads  one  to  conclude  that  it  is  not 
possible  to  describe  the  velocity  profile  in  the  outer  region  for  a  separated 
flow  by  the  universal  wake  function.  No  universal  backflow  function  appears 
to  exist. 

V. 2  Flow  Detachment  and  Upstream-Downstream  Intermittency 


It  is  well  established  that  separation  of  a  turbulent  boundary  layer  does 
not  occur  at  a  single  streamwise  location  but  is  spread  over  a  streamwise  region 
and  involves  a  spectrum  of  states.  Sandborn  and  Kline  (1961)  and  Sandborn  and 
Liu  (1968)  defined  the  limiting  points  of  the  region  as  the  "intermittent"  and 
the  "fully-developed"  separation  points.  The  former  indicates  the  onset  of  sepa¬ 
ration  by  the  appearance  of  intermittent  backflow  and  the  latter  signifies  the 
vanishing  of  the  mean  wall  shear  stress. 

Sandborn  and  Liu  (1968)  gave  correlations  between  H12  and  6^6  gg5  to 
demarcate  the  regions  of  intermittent  and  fully-developed  separation.  Figure  32 
gives  their  correlations  and  the  present  experimental  data  points.  According  to 
their  correlations,  the  present  data  show  intermittent  separation  to  occur  at 
130  inches.  The  value  of  YpU0  at  that  point  is  0.81  which  very  nearly  coincides 
with  the  value  obtained  by  Simpson  et  al.  (1977)  and  is  also  in  reasonable 
agreement  with  the  value  obtained  by  Sandborn  and  Liu.  By  interpolation  the 
fully-developed  separation  point  occurred  at  140  inches. 

At  the  recent  Project  SQUID  Colloquium  on  Flow  Separation  (Simpson,  1979), 
it  was  pointed  out  that  the  term  "separation"  must  mean  the  entire  process  of 
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Figure  32.  Sandborn's  H12  vs.  61 /5Q  995  separation  correlation: 

symbols.,  present  data;  solid  line,  intermittent  separation; 
broken  line,  fully-developed  separation. 


"departure"  or  "breakaway"  or  the  breakdown  of  boundary-layer  flow.  An 
abrupt  thickening  of  the  rotational  flow  region  next  to  a  wall  and  signi¬ 
ficant  values  of  the  normal -to-wall  velocity  component  must  accompany 
breakaway,  else  this  region  will  not  have  any  significant  interaction  with 
the  freestream  flow.  A  set  of  quantitative  definitions  were  proposed  and 
are  shown  on  Figure  33  along  with  old  definitions.  Figure  23  shows  the 
locations  of  incipient  detachment,  intermittent  transitory; detachment,  and 
transitory  detachment  for  the  present  flow  obtained  from  Figure  20.  In 
describing  a  quantitative  amount  of  backflow,  the  word  "detachment"  was 
preferred  over  "separation"  since  the  latter  term  refers  to  the  entire 
phenomenon.  Here  we  will  continue  to  use  the  time-honored  terminology,  but 
mention  the  new  terminology  for  the  sake  of  completeness. 

Downstream  of  intermittent  separation,  Simpson  et  al .  (1977)  showed 


the  existence  of  similarity  in  Yf 
pu  puo  vs.  y/M  where  y 


distributions  by  normalizing  and  plotting 
was  taken  as  the  value  near  the  wall  as  obtained 


from  a  figure  similar  to  Figure  20  and  M  was  the  distance  of  the  peak  in  the  u1 

distribution  from  the  wall.  The  present  data  also  exhibit  similarity,  particularly 

in  the  region  0.1  £  y/M  <  1.0,  with  it  improving  as  one  moves  downstream.  In 

fact  the  last  two  stations  at  156.4  inches  and  170.9  inches  show  the  similarity 

to  exist  all  across  the  boundary  layer,  including  the  backflow  region.  The 

similarity  in  the  backflow  region  improves  when  the  minimum  value  of  YpU  is 

used  instead  of  ydu  as  shown  in  Figure  34.  This  is  due  to  the  relatively  large 
puo 


uncertainty  in  y 


Simpson  et  al .  (1977)  curve-fitted  their  data  and  gave  an 


80 


Figure  34.  (y  -  y  )/0  -  Y  .  )  vs.  y/M;  y 

K  K  imn  K  min  K  min 

obtained  from  Figure  19;  M  is  distance  from  wall 
to  u'max”>  solid  line,  distribution  from  data  of 

Simpson  etal.  (1977). 


equation  for  the  distribution  in  the  region  0.1  £  y/M  <  1.0.  Figure  34 
shows  that  the  present  data  approximately  satisfy  the  equation.  Similar 
plots  drawn  with  M  being  taken  as  the  distance  from  the  wall  to  the  location 
where  peaks  were  observed  in  the  v'  and  -uv  distributions  show  as  good  or 
better  similarity,  such  as  in  Figure  35. 

Figures  36  (a)  and  (b)  show  results  forypv  or  the  fraction  of  time 

that  the  flow  is  away  from  the  wall.  Because  the  uncertainties  in  ypy  are 

relatively  large  near  the  wall,  y  was  used  in  the  normalized  results 

pvmin 

shown  in  Figure  37  for  the  region  downstream  of  intermittent  separation. 

Near  the  outer  edge  of  the  boundary  layer  the  intermittency  is  everywhere 
approximately  equal  to  one,  indicating  that  the  flow  is  always  directed  out¬ 
wards.  Near  the  wall,  the  intermittency  ypv  obtained  in  the  region  down¬ 
stream  of  intermittent  separation  is  higher  than  the  values  attained  upstream 
of  it,  which  can  be  attributed  to  the  flow  leaving  the  wall  as  a  consequence 
of  intermittent  separation.  As  in  the  case  of  ypu,  the  distributions  near 
the  wall  are  trough-shaped  in  the  region  downstream  of  intermittent  separation 
and  show  some  similarity. 

V.3  Turbulence  Correlations 

A.  Reynolds  Stresses  Correlations 

Figures  38  show  distributions  of  the  shear  stress  correlation  coefficient 
-uv/uV,  which  is  a  measure  of  the  extent  of  correlation  between  u  and  v 
fluctuations.  Table  5  gives  typical  uncertainty  values  for  the  correlation 
coefficients  presented  here  for  the  central  portion  of  the  boundary  layer. 

Near  the  outer  edge  the  values  are  larger  since  -uv,  u'  and  v'  approach  zero. 
Figure  38  (a)  also  shows  distributions  for  the  Schubauer  and  Klebanoff  (1951) 
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Figure  36(b). 


Fraction  of  time  that  the  flow  moves  away 
from  the  wall,  y_  vs.  y/6.  Lines  for  visual 
aid  only.  Note  displaced  ordinates. 
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for  the  86.5  inches  location.  As  recommended  by  Cebeci  and  Bradshaw  (1977),  a 
constant  value  of  0.08  is  used  for  £/<$  in  the  outer  region.  The  present  data  at 
86.5  inches  are  in  reasonable  agreement  with  these  results. 

Although  the  downstream  stations  exhibit  similarity  in  the  inner  layer,  they 
show  a  continuously  decreasing  mixing  length  in  the  outer  layer  as  one  moves  down¬ 
stream.  Further  downstream  in  the  intermittent  separation  region,  the  inner  layer 
similarity  gradually  disappears  and  the  mixing  length  in  the  outer  layer  continues 

to  decrease  with  no  region  of  constant  mixing  length.  In  the  separated  region, 

8  U 

Prandtl's  mixing  length  cannot  be  defined  in  the  backflow  region  where  gy  is 
negative.  The  distributions  for  the  forward  flow  region  are  shown  in  Figure  41  (d). 
They  indicate  large  values  of  the  mixing  length  closest  to  the  wall  where  it  can  be 
defined,  decreasing  continuously  as  one  moves  farther  away  from  the  wall.  There  is 
also  some  indication  of  the  profiles  achieving  similarity. 

Figure  42  show  the  eddy  viscosity  profiles  in  the  various  regions.  As  in  the 
case  of  the  mixing  length,  a  few  sets  of  data  from  earlier  investigations  are  also 
plotted  for  comparison.  In  general,  the  same  comnents  made  about  the  mixing  length 
profiles  are  applicable  to  these  profiles  also.  The  present  data  in  Figure  42  (a) 
show  good  agreement  with  Klebanoff's  (1955)  data  in  the  zero  pressure  gradient 
region.  The  data  in  Figure  42  (b)  show  good  agreement  with  Bradshaw's  data  in  the 
adverse  pressure  gradient  region  in  the  inner  layer.  A  prediction  using  Cebeci  and 
Smith's  model  in  the  relation 

Ve  ■  *2  |  (23) 
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TABLE  6:  Typical  Uncertainties  for  turbulence  corrections 


Location 

Quantity 

Type  of  Data 

Absolute  Value 

Absolute  Value  of 

X 

Y 

tne  uncertainty 

percentage  un- 

inches 

inches 

at  tne  particular 

certainty  for 

x-y  location 

the  complete 
data  set 

112.375 

0.5 

-uv 

L.D.V. 

0.462 

0.031 

6.6 

111.25 

0.6 

v-V 

x-wire 

0.399 

0.076 

19 

112.375 

0.5 

-uv 

L.D.V. 

0.204 

0.014 

6.8 

111.25 

0-6 

(vW\) 

x-wire 

0.174 

0.034 

19.5 

112.375 

0.5 

tan* -2  vlv  ^  j 

L.D.V. 

40.9° 

1 .8® 

4.4 

111.25 

0.6 

u.'  -v' 

x-wire 

35,4" 

6.5" 

18.4 

131.875 

1.0 

,  ,2  ,2.  3U 

-<u  -v  >  ft 

Smoothed 

X-wire  4 

0.389 

0.132 

35 

L.D.V. 

_ 2J _ 

118.5 

2.0 

2  2 
-V  ) 

Smoothed 

-0.411 

0-314 

Varies  widely 

131.875 

2.5 

!>-) 

L.D.V. 

-0.742 

0.041 

131.875 

1.5 

t. 

L.D.V. 

0.052 

0.0056 

12 

86.5 

x-wire 

0.064 

0.013 

20 

131.875 

1.5 

L.D.V. 

0.006184 

0.00074 

15 

86.5 

0.354 

_ 

x-wi re 

0.0133 

0.0027 

20 

Table  6:  Flow  conditions  for  the  present  and  previous  Investigations. 


streamwise  location 

Data  of  other  investigators 

Parameter 

(in  inches)  for  the 
present  data 

streamwise 
location  (in 

ft)  for  Bradshaw 

Klebanoff ‘s 

data 

East  i 
Sawyer, 
flow  4 

East  6 
Sawyer, 
flow  1 

86.5 

105 

117.6 

17.5  22.5  a*-0.15  a*-0.255 

8  (Inches) 

0.153 

0.284 

0.458 

0.338 

0.129 

Re 

5205 

8617 

11988 

18750  41850  22900  38800 

H 


12 


1.4 


1 .418  1.625  2.024 


1.35  1.6 


1.54 


1.344 


1.31 


****** 1 


*3 


strong  adverse  gradient  boundary  layer.  These  two  sets  of  measurements 
compare  reasonably  well,  considering  the  fact  that  the  adverse  pressure 
gradient  distributions  are  different.  Table  6  shows  a  comparison  of 
some  parameters  for  the  two  flows. 

Figure  38  (b)  shows  distributions  in  the  vicinity  of  the  beginning 
of  intermittent  backflow.  Unlike  the  distributions  far  upstream  shown  in 
Figure  38  (a)  or  those  observed  in  zero  pressure  gradient  boundary  layers, 
the  distributions  in  this  region  do  not  exhibit  a  constant  value  over  a 
large  part  of  the  outer  layer.  However,  the  distributions  for  some  of  the 
stations  do  indicate  a  small  region  with  a  nearly  constant  value  as  low  as  0.2 
to  0.3.  As  one  moves  downstream,  the  peaks  for  the  distributions  seem  to 
gradually  move  towards  the  outer  edge  of  the  boundary  layer.  Similar  features 
such  as  correlation  coefficients  as  low  as  0.3  with  the  peaks  occuring  near  the 
outer  edge  of  the  boundary  layer  were  observed  by  Spangenberg  et  al.  (1967)  in 
their  experiments  on  an  adverse  pressure  gradient  flow  approaching  separation. 
Not  much  significance  can  be  attached  to  the  dips  in  the  distributions  observed 
near  the  wall  except  to  hint  that  they  might  be  a  consequence  of  the  peaks  in 
the  production  curves  occurring  near  the  wall.  Figure  38  (c)  indicates  that 
the  profiles  for  the  separated  region  seem  to  exhibit  some  similarity.  These 
distributions  compare  fairly  well  in  the  outer  region  with  the  results  of 
Wygnanski  and  Fielder  (1970)  for  a  mixing  layer.  Figure  39  gives  the  distri¬ 
butions  of  another  type  of  correlation  coefficient,  a^  =  -u\^u'^  +  v’^),  similar 

to  the  one  used  by  Bradshaw  et  al .  (1967)  for  converting  the  turbulent  kinetic 

2  12  2 

energy  equation  into  an  equation  for  shear  stress.  Using  w'c  =  j  (u*  +  v’  ),  it 

is  possible  to  relate  to  the  more  commonly  used  Bradshaw's  constant  'a' 

defined  as  a  -  -uv  /  (u1  +  v'  +  w'  )  by  the  relation  a  =  2/3af. 
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Figure  38(a).  Shear  stress  correlation  coefficient  profiles 
from  smoothed  laser  and  hot-wire  data.  Figures 
16,  -uv/u'v'  vs.  y/6  :  upstream  of  separation. 
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Figure  38(b).  Shear  stress  correlation  coefficient  profiles 
from  smoothed  laser  and  hot-wire  data.  Figures 
16,  -uv/u'v'  vs.  y/6:  in  the  vicinity  of  the 
beginning  of  intermittent  backflow. 
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Figure  38(c).  Shear  stress  correlation  coefficient  profiles  from 
smoothed  laser  and  hot-wire  data.  Figures  16, 
-uv/uV  vs.  y/6:  separated  region. 


I  Figure  39(a).  Shear  stress  correlation  coefficient  profiles 

from  smoothed  laser  and  hot-wire  data,  Figures 

116,  -uv/  (u ' 2  +  v'^)  vs.  y/5:  upstream  of 

separation. 

I 
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Figure  39  (a)  also  contains  the  data  of  East  and  Sawyer  (1979)  for 
favorable  and  adverse  pressure  gradient  flows.  The  flow  conditions  for  those 
cases  are  given  in  Table  6.  Considering  the  wide  variations  in  the  flow  condi¬ 
tions  and  the  uncertainties  in  the  measurements,  the  agreement  seems  to  be 
reasonable,  particularly  for  the  adverse  pressure  gradient  case.  The  variation 
in  the  behavior  of  the  distributions  as  one  moves  downstream  is  similar  to  that 
for  the  shear  correlation  coefficient  -uv/u'v1,  with  an  increasingly  reduced 
flat  region  and  a  reduction  -n  the  value  of  a-j  to  as  low  as  0.1  for  the 
separated  region. 

_  1  "1 

Another  quantity  which  can  be  derived  from  u1,  v'  and  -uv  is  0  =  j  tan 
-2uv  .  which  gives  the  angle  of  inclination  of  the  principal  axis  to  the 
u,2-v'2 

flow  direction.  This  has  been  plotted  in  Figures  40.  Sandborn  and  Slogar  (1955) 
observed  that  e  remains  almost  independent  of  x  and  y  except  for  a  small  part  of 
the  inner  layer.  They  also  noticed  that  in  the  inner  layer  the  angle  6  decreases 
rapidly  as  the  wall  is  approached,  thus  tending  to  align  the  axis  of  the  principal 
stress  with  the  flow  direction.  Considering  the  uncertainties  in  6 ,  particularly 
near  the  outer  edge  where  all  the  quantities  u',  v1  and  -u7-*>0,  the  present  data 
seem  to  indicate  those  same  trends,  at  least  for  the  stations  downstream  of  86.5 
inches  and  up  to  the  beginning  of  the  intermittent  separation.  In  the  inter¬ 
mittent  separation  region,  only  some  of  the  stations  indicate  a  weak  dependence  with 
respect  to  y  in  the  outer  layer.  In  the  fully-separated  region  there  is  an 
indication  of  the  profiles  tending  to  become  similar.  The  angle  in  the  flat 
region  decreases  from  approximately  18°  at  86.5  inches  to  12°  in  the  separated 
region. 
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Figure  40(a).  Angle  of  principal  axis  of  stress  to  the  ? 

flow  direction,  0  =  1/2  tan(-2uv/(u'2  -  v'^)) 
vs.  y/5.  Solid  line  is  mean  line  from  Sandborn  and 
Slogar  (1955). 


Figure  40(b).  Angle  of  principal  axis  of  stress  to  the 

flow  direction,  0  =  1/2  tan(-2uv/u'^  -  v'^)) 
vs.  y/5.  Solid  line  is  mean  line  from  Sandborn 
and  Slogar  (1955). 
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Figure  40(c).  Angle  of  principal  axis  to  the  flow 

direction,  0  =  1/2  tan(-2uv/(u'2  -  v’2)) 
vs.  y/6.  Solid  line  is  mean  line  from 
Sandborn  and  Slogar  (1955). 


B.  Eddy  viscosity  and  Prandtl  mixing  length  distributions 
The  Prandtl  mixing  length 


l  -UV  Bill  /  3U  \ 
6  '  6  3y  / 


and  the  eddy  viscosity 


(20) 


(21) 


were  calculated  from  measured  Reynolds  shearing  stress  and  calculated  velocity 
gradient  distributions.  Figure  41  (a)  shows  the  mixing  length  results  for  the 
region  up  to  the  throat  of  the  test  section  where  the  pressure  gradient  is 
either  favorable  or  approximately  zero.  The  data  of  Klebanoff  (1955)  for  a 
zero  pr  ssure  gradient  boundary  layer  and  that  of  East  and  Sawyer  (1979)  for 
zero  and  favorable  pressure  gradient  boundary  layers  are  also  presented  for 
comparison.  The  present  data  at  64.2  inches  show  good  agreement  within  the 
limits  of  uncertainty  with  the  zero  pressure  gradient  data  of  the  earlier 
investigators.  The  data  at  31.25  inches  show  agreement  only  in  the  inner  layer 
with  the  favorable  pressure  gradient  data  of  East  and  Sawyer.  One  possible 
reason  for  this  might  be  the  close  proximity  of  that  station  to  the  entrance 
region  of  the  test  section. 

Figure  41  (b)  covers  the  adverse  pressure  gradient  region  of  the  flow  up 
to  the  start  of  the  intermittent  separation.  The  data  of  Bradshaw  (1967)  for 
adverse  pressure  gradient  equilibrium  boundary  layers  and  East  and  Sawyer  (1979) 
are  presented  for  comparison.  Also  shown  in  Cebeci  and  Smith's  (1974)  extension 
of  van  Driest’s  mixing  length  model  for  the  inner  layer 
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Figure  41(d).  Mixing  length  distributions,  1/ 5  vs.  y/6 
downstream  of  separation. 
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Table  7:  The  ratio  of  UT<Sgg  for  successive  locations 


Streamwise 
location  in 
inches 

Ut699 

ft.2/sec. 

Ratio  of  U  6__ 
t  yy 

for  successive 
stations 

86.5 

0.255 

1.066 

1.058 

n  OfiQ 

105.3 

0.272 

112.4 

0.288 

• 

118.5 

0.262 

is  in  reasonable  agreement  with  the  inner  layer  data  at  86.5  inches. 

At  first  it  is  a  little  surprising  that  there  is  similarity  in  the 
inner  layer  mixing  length  distributions  and  similarity  in  the  inner  layer 
eddy  viscosity  distributions  near  separation  when  6  is  used  for  scaling  y. 

However,  as  shown  in  Table  7  the  ratio  of  UTS  at  successive  stations  is 
near  unity  in  this  region,  so  y+/(y/6)  is  the  same  for  successive  stations 
and  the  profiles  near  the  wall  are  similar  with  respect  to  y+  as  well.  In 
the  intermittent  separation  region,  the  inner  layer  similarity  disappears 
and  the  eddy  viscosity  decreases  with  respect  to  x  in  the  outer  layer.  In 
the  separated  region,  vg  can  be  defined  everywhere  except  where  3U/8y  =  0. 

Eddy  viscosity  profiles  also  show  some  similarity  in  the  outer  layer  as  well 
as  near  the  wall  in  the  separated  region. 

For  both  mixing  length  and  eddy  viscosity,  the  data  in  the  vicinity  of 
separation  indicate  much  lower  values  in  the  outer  region  than  for  attached 
boundary  layers.  As  shown  below  in  section  V.5,  normal  stresses  effects  can 
be  used  to  explain  this  behavior. 

C.  Skewness  and  flatness  factor  distributions 

Some  measurements  of  skewness  and  flatness  factors  of  the  u  and  v  fluctua¬ 
tions  have  been  done  in  zero  pressure  gradient  boundary  layers  and  in  channel 
flows  by  Dumas  (1966),  Zaric  (1972),  Kreplin  (1973),  Antonia  (1973)  and  lleda 
and  Hinze  (1975).  Only  Sandborn  (1959)  is  known  to  have  made  measurements  of 
the  flatness  factor  in  an  adverse  pressure  gradient  boundary  layer  flow  in  the 
vicinity  of  separation. 

Figures  43  (a)  and  44  (a)  show  a  comparison  of  the  present  laser  anemometer 
data  for  Fu  and  Fy  with  the  zero  pressure  gradient  boundary  layer  data  of  Antonia 
(1973).  The  good  agreement  observed  between  the  two  sets  of  data  in  the  logarithmic 


•  Sandborn;  Turbulent  bohoary  layer  in  the 
vicinity  OF  separation;  R^“  5687 
a  Antonia;  Turbulent  bcxmary  layer;  Re*  31000 

▼  Dumas;  Turbulent  boundary  layer;  Rf»  32500 
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Figure  43(a).  Flatness  factor  F  profiles:  upstream 
results.  1 
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Figure  44(b).  Flatness  factor  F  profiles:  downstream 
of  separation.  v 
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region  and  the  outer  region  indicates  that  the  pressure  gradient  does  not 


have  much  effect  on  Fu  and  Fy  in  those  regions.  Comparison  with  Figures  43  (b) 
and  44  (b)  for  the  flow  downstream  of  separation  indicates  that  separation  also 
does  not  have  much  effect  on  Fy  and  Fy  over  the  shear  layer. 

However,  when  plotted  against  y+  in  Figure  45,  the  data  for  Fu  upstream  of 
separation  indicate  an  apparent  effect  of  pressure  gradient  in  the  region  close 
to  the  wall,  mainly  in  the  buffer  layer.  In  the  viscous  sublayer  for  both  zero 
and  adverse  pressure  gradient  flows,  the  flatness  factor  attains  values  much 
higher  than  the  value  for  a  gaussian  probability  distribution,  which  is  equal  to 
3.  This  is  possible  because  the  inrush  phase  of  the  bursting  cycle  which  brings 
in  high  velocity  fluid  from  the  outer  region  results  in  large  amplitude  positive 
u  fluctuations  and  consequently  produces  a  large  skirt  in  the  velocity  probability 
distribution.  Similarly,  near  the  outer  edge  of  the  boundary  layer,  intermittent 
large  amplitude  negative  u  fluctuations  occur  as  a  result  of  the  large  eddies 
driving  the  fluid  from  the  low  velocity  regions  outwards,  which  tends  to  increase 
the  flatness  factor.  In  the  buffer  layer  near  a  y+  of  13,  the  zero  pressure 
gradient  flows  of  Ueda  and  Hinze,  Zaric,  and  Kreplin  all  show  a  dip  in  the  Fu 
flatness  factor  distributions  and  a  change  in  sign  in  the  skewness  factor  Sy 
distributions  for  u  as  shown  in  Figure  46.  Ueda  and  Hinze  have  remarked  that  this 
location  is  where  u'  attains  the  maximum  value.  The  present  data  neither  show 
any  such  predominant  dip  in  Fy  nor  sign  change  of  Su  in  the  buffer  layer. 
Sandborn's  (1959)  data  for  Fu  in  an  adverse  pressure  gradient  boundary  layer 
flow  show  a  behavior  similar  to  the  present  data.  The  present  data  for  Fu  and 
Su  also  show  reasonable  agreement  with  those  of  Dumas  (1966),  but  the  significance 
of  this  is  clouded  since  the  pressure-gradient-flow  conditions  were  not  mentioned 
in  his  paper. 
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The  present  data  for  as  shown  in  Figure  47  (a)  indicate  a  change 
in  sign  at  a  location  farther  away  from  the  wall  (y/6  ~  0.4).  This  location 
corresponds  to  the  region  where  the  Reynolds  shear  stress  and  the  turbulent 
intensities  reach  their  maximum  values.  The  intense  momentum  exchange  in 
this  region  results  in  the  lack  of  occasionally  very  high  or  very  low  fluctua¬ 
tions  and  as  a  consequence  the  probability  distribution  does  not  have  much 
skewness.  As  one  moves  closer  to  the  wall,  the  intermittent  large  amplitude 
positive  u  fluctuations  tend  to  make  the  probability  distributions  more  posi¬ 
tively  skewed  (Eckelmann,  1974)  and  vice-versa  when  one  moves  away  from  the 
wall . 

The  location  corresponding  to  zero  skewness  for  u  occurs  very  close  to 
the  wall  in  zero-pressure-gradient  flows  because  the  Reynolds  shear  stress 
attains  a  maximum  value  in  that  region.  Furthermore,  the  intense  mixing  in 
that  region  surpresses  large  amplitude  u  fluctuations,  thus  removing  the  skirt 
in  the  positively  skewed  velocity  probability  distribution  and  changing  it  to 
a  more  nearly  top-hat  shape  with  a  low  flatness  factor.  The  same  does  not 
happen  in  adverse  pressure  gradient  flows  in  the  region  of  maximum  shear  because 
the  probability  distribution  in  that  region  is  more  nearly  gaussian  with  only  a 
slight  skewness  and  with  no  significant  large  amplitude  fluctuations  to  be 
suppressed. 

Downstream  of  separation  the  skewness  Su  is  reduced  to  negative  values  in 
the  backflow  region  as  shown  in  Figure  47  (b).  A  maximum  is  observed  in  the 
vicinity  of  the  minimum  mean  velocity.  As  shown  in  Figure  43  (b),  Fu  also  has 
a  local  maximum  near  this  location.  The  second  zero-skewness  point  is  slightly 
closer  to  the  wall  than  the  location  of  the  maximum  shear  stress. 

The  flatness  factor  distributions  for  v  in  Figures  44  (a)  and  (b)  show  a 
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Figure  47(b).  Skewness  factor  S  vs.  y/6  profiles  laser 
anemometer  data:  results  downstream  of 
separation. 


trend  similar  to  that  of  u,  the  only  difference  being  the  reduced  width 
of  the  flat  region.  Figures  48  and  49  show  that  there  is  a  significant 
variation  of  Sy  along  the  flow.  Only  downstream  of  112  inches  is  there 
profile  similarity  in  the  outer  region.  Sy  shown  in  Figure  49  exhibit  a 
shape  approximately  opposite  in  sign  to  that  of  5y,  with  a  large  positive 
skewness  factor  near  the  outer  edge  of  the  boundary  layer,  gradually 
decreasing  to  negative  values  towards  the  wall.  This  results  in  the  appear¬ 
ance  of  two  zero-skewness  points  in  the  distributions  of  Sy  both  upstream 
and  downstream  of  separation.  The  zero-skewness  point  which  is  farther  from 
the  wall  occurs  in  the  region  of  maximum  shear  both  upstream  and  downstream 
of  separation,  which  indicates  that  the  backflow  has  no  influence  on  the 
location  of  this  point  as  in  the  case  of  Su-  Downstream  of  separation  the 
flatness  and  skewness  factors  away  from  the  wall  are  in  qualitative  agreement 
with  those  of  Wygnanski  and  Fiedler  (1970)  for  a  plane  mixing  layer.  This  is 
not  surprising  since  the  mean  velocity  profiles  resemble  those  in  mixing 
layers. 

D.  Diffusion  of  turbulence  kinetic  energy 

_  ~2~ 

The  diffusion  term  9/9y  (pv/p  +  1/2  q  v)  of  the  turbulence  kinetic  energy 
equation  is  known  to  become  more  important  as  a  turbulent  boundary  layer  approaches 
separation  (Bradshaw,  1967b;  Simpson  and  Collins,  1978).  The  term  pv/p  which 
represents  the  diffusion  flux  due  to  pressure  forces  cannot  be  measured  directly 
using  available  techniques.  Normally,  it  is  estimated  by  the  difference  of  other 
measureable  terms  in  the  turbulence  kinetic  energy  equation,  although  experi¬ 
mental  uncertainties  make  the  results  quite  uncertain.  Here  the  turbulence 
kinetic  energy  diffusion  flux 
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Figure  98.  Skewness  factor  Sy  vs.  y/«  profile  from  cross  hot-wire  anejocter  data 
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Figure  49(a).  Skewness  factor  S  vs.  y/6  from  laser 
anemometer  data:  vupstream  results. 
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Figure  49(b).  Skewness  factor  S  vs.  y/6  from  laser 
anemometer  data:  results  downstream  of 
separation. 
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“5“  ~ 3 

was  estimated  using  u  v  and  v  cross-wire  anemometer  measurements  and  the 

~2~  ~T~  ~ 3 

approximation  proposed  by  Bradshaw  (1967b),  w  v  =  (u  v  +  v  )/2. 

Figures  50  (a)  and  (b)  show  the  present  results.  The  flux  of  turbulence 
kinetic  energy  is  positive  in  the  regions  where  data  have  been  plotted,  indi¬ 
cating  that  the  flux  is  directed  away  from  the  wall.  For  locations  downstream 
of  117.6  inches  the  data  are  limited  only  to  the  region  near  the  outer  edge  of 
the  boundary  layer.  Nearer  the  wall  at  these  locations  the  flux  is  expected  to 
be  negative,  since  most  of  the  turbulence  energy  production  is  in  the  middle  of 
the  boundary  layer  and  previous  strong  adverse  pressure  gradient  data  (East  and 
Sawyer,  1979)  have  this  behavior. 

East  and  Sawyer  proposed  a  gradient  model  based  on  a  mixing  length 
formulation 

3/2 

=  0.4  A  ^  (7)  (25) 

They  obtained  experimental  data  for  seven  equilibrium  turbulent  boundary  layers 
with  U  ~  xR  and  R  approximately  equal  to  0.4,  0.2,  0,  -0.2,  -0.4,  -0.6,  and  -0.8. 
The  above  model  agreed  with  those  data  satisfactorily  in  the  outer  half  of  the 
boundary  layer  in  all  cases.  Agreement  in  the  inner  regions  improved  for  increas¬ 
ingly  adverse  pressure  gradients.  Using  the  mixing  length  and  turbulence  kinetic 
energy  distributions  obtained  from  the  present  equilibrium  experiments,  the  results 
from  this  model  are  shown  in  Figures  50.  Agreement  in  the  outer  region  is  within 
the  uncertainty  of  the  measurements.  In  the  inner  region  only  the  general  shape 
of  the  predictions  agree  with  measurements. 

It  can  be  observed  from  Figures  51  (a)  and  (b)  that  the  diffusion  is  small 
at  the  upstream  stations,  becoming  appreciable  downstream  from  117  inches.  Farther 
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Figure  50(a).  Turbulent  kinetic  energy  diffusion  flux 

vs.  y/6  -  hot-wire  anemometer  results  com¬ 
pared  to  East  and  Sawyer's  gradient  diffusion 
model:  upstream  of  separation. 
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Figure  51(a).  Turbulent  kinetic  energy  diffusion  vs.  y/S 
downstream  of  separation. 
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Figure  51(b).  Turbulent  kinetic  energy  diffusion  vs 
y/6:  downstream  of  separation. 


downstream  as  separation  is  approached,  the  diffusion  increases  continuously. 

It  is  interesting  to  note  that  such  large  negative  diffusion  rates  occur  on 
the  low  velocity  side  of  mixing  layers  also.  This  can  be  seen  in  Fig.  51  (b) 
which  has  the  data  of  Wygnanski  and  Fiedler  (1970)  plotted  for  comparison  with 
the  present  data  at  x  =  156  3/8  inches.  The  maximum  velocity  in  the  mixing 
layer  Um  and  the  total  shear  layer  thickness  2y^  were  used  for  nondimension- 
alizing  those  data.  This  similarity  in  behavior  with  the  mixing  layer  suggests 
that  the  diffusion,  which  is  responsible  for  the  lateral  spread  of  mixing  layers, 
is  also  responsible  for  the  rapid  growth  of  separated  boundary  layers.  The 
large  gain  of  energy  by  diffusion  in  the  outer  region  and  the  associated  increase 
in  entrainment  of  the  nonturbulent  fluid  seems  to  be  responsible  for  the  mainten¬ 
ance  of  the  large  eddies  and  the  large  growth  rates  of  separated  boundary  layers. 

The  increase  in  entrainment  rate  of  free-stream  fluid  as  separation  is 
approached  is  demonstrated  in  Figure  52  in  terms  of  the  entrainment  velocity,  Vp, 
obtained  from  mean  velocity  measurements  using  the  relationship 

V, 


=  dx  _U»  (60.995  " 


) 


(26) 


Upstream  of  separation  these  results  are  in  good  agreement  with  Bradshaw's  (1967) 
correlation 
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for  boundary  layer  and  mixing  layers.  Downstream  of  separation  there  is  poor 
agreement,  in  contrast  to  the  good  agreement  obtained  by  Simpson  et  al.  (1977) 
for  their  separating  flow.  This  might  be  because  of  some  three-dimensionality 
which  seems  to  exist  in  that  region. 


Figure  53  shows  the  distribution  of  the  diffusion  function  G/(t  /p)2, 
which  was  defined  by  Bradshaw  (1967a)  to  relate  the  turbulence  kinetic  energy 
diffusion  to  the  turbulent  shear  stress. 

s  -  (py/p  t.  h  ,q2»J.  (28) 

(r/p) (t  /p) 2 
'  max  ' 


The  diffusion  function  used  by  Bradshaw  and  that  computed  from  the  data  of 
East  and  Sawyer  (1979)  are  also  shown.  Although  there  are  large  differences 
up  to  half  of  the  boundary  layer  thickness,  the  present  data  blend  in  with  their 
data  in  the  outer  region.  When  compared  with  Bradshaw’s  diffusion  function,  the 
differences  are  larger  and  there  is  no  region  of  agreement  of  all.  The  diffusion 
function  given  by  Bradshaw  was  derived  from  the  zero  pressure  gradient  boundary 
layer  data  of  Klebanoff  (1955).  These  results  indicate  that  the  diffusion  function 
is  dependent  upon  pressure  gradient  conditions. 

V.4  Momentum  and  Turbulence  Energy  Balances 

In  order  to  further  understand  the  effect  of  separation  on  the  transport  of 
momentum  and  turbulence  kinetic  energy,  terms  of  the  governing  equations  were  obtained 
using  the  measured  quantities  described  above.  The  x-direction  and  y-direction 
momentum  equations  are,  respectively 


u  9U  +  v  9U  _-l  3P  +  3(-uv)  .  3? 

u  3x  v  9y  ~  p  9x  9y  9x 


(29) 


9V  v  9V  _-l  9P  +  9(-uv)  _  _9v_ 
9x  *  3y  "  p  9y  9x  ”  9y 


(30) 


For  each  equation  the  terms  on  the  left  side  are  inertia  or  convective  terms  while 
the  terms  on  the  right  side  describe  the  pressure  gradient,  the  shearing  stress 
gradient,  and  the  normal  stress  gradient,  respectively.  The  turbulence  energy 
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Figure  53.  Bradshaw's  (1967)  diffusion  parameter 
vs.  y/6  g95  for  several  sets  of  data. 
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equation  is 


U  3q2  V  3q2  _ 

2  ax  2  ay 


-a 


V2)  f  *  £ 


(31) 


The  terms  on  the  left  side  are  advection  terms  while  the  terms  on  the  right 
side  describe  turbulent  diffusion,  turbulent  shear  stress  production,  normal 
stresses  production,  and  dissipation,  respectively.  Dissipation  was  not 
measured.  In  all  three  equations  the  viscous  terms  have  been  neglected  since 
they  are  much  smaller  than  the  other  terms. 

An  estimate  of  uncertainties  of  all  the  significant  terms  for  a  few 
typical  points  across  the  boundary  layer  are  given  in  Table  8  for  118.5  inches 
and  131.875  inches.  Very  near  the  wall  the  uncertainties  are  high,  but  beyond 
y/<$  x  0.02,  the  uncertainties  of  most  of  the  dominant  terms  are  less  than  30% 
40%  at  many  points.  In  general,  the  terms  involving  derivatives  with  respect 
to  y  have  less  uncertainty  as  compared  to  those  involving  derivatives  with 
respect  to  x,  since  the  latter  terms  are  much  smaller  and  were  computed  from 
data  acquired  on  different  days.  Hence  each  data  point  used  to  determine  x 
derivatives  corresponded  to  slightly  different  experimental  conditions. 

An  exception  to  this  is  the  inertia  terms  of  the  x-di recti  on  momentum 
equation.  In  this  case  the  two-dimensional  continuity  equation  can  be  used  to 
obtain  a  single  term  involving  only  a  y  derivative  of  a  given  velocity  profile 


U  3U  +  V  =  -U2  3  Mu) 

u  9x  v  3y  u  3y 


I 


This  expression  was  used  only  when  U  was  much  larger  than  V,  since  the  uncer¬ 
tainty  in  V/U  becomes  large  as  U  approaches  zero.  The  relative  uncertainty  in 
this  term  is  large  in  the  outer  region  because  y-direction  gradients  are  small 
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Estimate  of  Uncertainties  for  the  terms  of  the  Momentum  and  Energy  eouatlons. 
Table  8a.  Terms  Involving  derivatives  with  respect  to  x  at  *  *  118.5  Inches. 


Estimate  of  Uncertainties  for  the  terms  of  the  Momentum  and  Energy  equations. 
Table  8c.  Terms  involving  derivatives  with  respect  to  x  at  x  1  131.875  inches. 
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Estimate  of  Uncertainties  for  the  terms  of  the  Momentum  and  Energy  eQuatlons. 
Table  8d.  Terms  Involving  derivatives  with  hesoect  to  y  at  x  *  131.875  Inches. 
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On  the  whole,  even  though  the  uncertainties  are  large  it  is  still  possible 
to  arrive  at  certain  conclusions  regarding  the  relative  importance  of  the 
various  terms  in  the  momentum  and  turbulence  energy  equations  as  the  boundary 
layer  passes  through  separation. 

Although  the  momentum  and  energy  balances  were  examined  at  a  number  of 
stations,  the  results  are  presented  here  for  three  representative  stations  only. 
They  correspond  to  a  location  upstream  of  separation  (118.5  inches),  a  location 
in  the  intermittent  separation  region  (131  7/8  inches),  and  one  in  the  fully- 
separated  region  (156  3/8  inches).  Figures  54  and  55  show  the  distributions  of 
the  various  non-dimensional  terms  of  the  momentum  equations  and  Figure  57  repre¬ 
sent  the  terms  of  the  energy  equation.  The  locations  of  the  maximum  shear  stress 
-uv  and  the  maximum  (u1  +  v'  )  are  shown  on  all  the  plots. 

Figure  55  indicates  that  the  only  important  terms  in  the  equation  for 
momentum  transport  in  the  y-di recti  on  are  the  pressure  gradient  and  the  normal 
stress  terms.  This  is  true  both  upstream  and  downstream  of  separation  and  leads 
to  the  following  simplification  of  eqn.  (30): 

-lit  -  av'2 
p  3y  9y 


2 

Upon  integration  it  becomes  P  (x,y)  =  -  pv'  .  Differentiating  this  equation 

with  respect  to  x  produces 


-1  9P  =  -1  3P°°  ,  3v'2 
p  8x  ’  p  3x  3y 


(32) 


The  pressure  gradient  evaluated  using  eqn.  (32)  is  also  plotted  in  Figure  54. 

3P 

A  first  look  at  Figure  54  indicates  large  discrepancies  between  ^  computed  using 
eqns.  (29)  and  (32).  However,  in  view  of  the  uncertainties  of  3P/9x  derived  from 
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eqn.  (29),  the  results  are  in  agreements  within  these  uncertainties. 

A  comparison  of  Figure  54  indicates  that  in  the  separated  region  the 
convective  terms  become  unimportant  in  the  inner  layer.  The  momentum  trans¬ 
fer  due  to  shear  mainly  balances  the  x-direction  pressure  gradient.  In  the 
outer  region  in  addition  to  the  important  convective  terms,  the  normal 
stresses  term  becomes  important  as  separation  is  approached,  as  has  already 
been  shown  by  Simpson  et  al .  (1977).  The  normal  stresses  play  an  important 
role  in  the  vicinity  of  the  maximum  shear  stress.  At  118  inches,  the  normal 
stresses  term  is  still  quite  small.  The  momentum  balance  at  112  inches  shows 
that  the  normal  stresses  term  is  more  important.  Its  importance  increases 
progressively  downstream  as  can  be  seem  from  Figures  54  (b)  and  (c),  which 
show  that  this  term  contributes  up  to  half  of  the  momentum  transport  in  the 
outer  region.  This  is  shown  more  clearly  in  Figures  56  (a)  and  (b)  by  the 
distributions  of  the  ratio  of  the  normal  stresses  term  to  the  shear  stress  term. 
However,  due  to  uncertainties  in  the  gradients  the  uncertainty  of  these  results 
in  the  outer  region  is  large,  as  shown  in  Table  5.  Thus  the  inner  layer  in  the 
separated  region  could  be  modeled  by  neglecting  the  convective  terms  while  in 
the  outer  layer  the  additional  effect  of  the  normal  stresses  must  be  included. 

Figures  57  show  the  importance  of  the  normal  stresses  turbulence  energy 
production  from  just  upstream  of  intermittent  separation  to  far  downstream.  The 
results  for  the  Bradshaw  (1967)  flow  are  in  qualitative  agreement  with  the  data 
shown  in  Figure  57.  (a).  Figures  58  (a)  and  (b)  show  the  ratio  of  the  normal 
stresses  production  to  the  shear  production  for  the  several  locations  in  the 
vicinity  of  separation.  As  indicated  by  the  present  data  and  the  data  of 
Simpson  et  al .  (1977)  and  Schubauer  and  Klebanoff  (1950),  the  normal  stresses 
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Figure  56(a),  Ratio  of  norma)  stresses  to  shear  stress 
terms  in  the  momentum  equation:  upstream 
of  separation. 
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effect  becomes  increasingly  important  as  separation  is  approached.  In  fact 
both  sets  of  SMU  data  show  good  agreement  in  the  corresponding  regions  of 
development,  with  a  near  doubling  of  the  ratio  in  the  intermittent  separation 
region.  The  present  data  in  that  region  indicate  the  presence  of  a  hump  in 
the  distributions  near  y/6  of  0.05  to  0.1,  which  becomes  more  significant  as 
separation  is  approached.  This  is  a  result  of  the  mean  velocity  profiles 
becoming  inflexional  in  nature,  which  produces  a  reduced  9U/3y  in  that  region. 

In  fact  these  humps  increase  rapidly  along  the  flow  until  3U/9y  attains  a  zero 
value  for  each  profile  in  the  backflow  region  where  the  velocity  reaches  a 
minimum  value.  The  earlier  data  of  Simpson  et  al .  (1977)  at  124.3  inches  also 
suggest  the  presence  of  a  hump.  In  the  backflow  region  the  two  types  of  pro¬ 
duction  oppose  each  other  as  shown  in  Figure  58  (b),  but  they  aid  one  another 
in  the  forward  flow  region.  The  distributions  in  the  outer  layer  tend  toward 
similarity  and  the  ratio  seems  to  be  almost  a  constant  of  0.6  for  0.2  <  y/6  <_  0.7. 

As  far  as  shear  production  alone  is  concerned,  the  present  data  in  the 
region  upstream  of  separation  is  in  agreement  with  those  of  Spangenberg  et  al, 
(1967)  and  others  who  observed  two  peaks  in  distributions  for  boundary  layers 
subjected  to  large  adverse  pressure  gradients.  The  present  data  indicate  that  as 
separation  is  approached,  the  peak  near  the  wall  becomes  weaker  until  it  vanishes 
in  the  region  of  fully-developed  separation.  In  the  backflow  zone  of  the  separ¬ 
ated  region  there  is  no  shear  production  as  indicated  by  Figure  57  (c)  and 
advection  is  also  insignificant.  Hence  the  only  mode  by  which  turbulence  energy 
can  reach  the  backflow  zone  is  by  turbulent  diffusion.  This  conclusion  is 
consistent  with  the  results  discussed  in  section  IV. 3. D  above:  diffusion  plays 
a  major  role  in  transporting  the  turbulent  kinetic  energy  in  separated  flows  from 
the  middle  part  of  the  layer,  where  it  is  mainly  produced,  to  the  outer  region  and 
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the  region  near  the  wall.  The  absence  of  production  near  the  wall  in 
separated  flow  also  leads  one  to  conclude  that  the  backflow  near  the  wall 
is  controlled  by  the  large-scaled  outer  region  flow,  rather  than  by  some 
wall -shear-stress-related  "law  of  the  wall". 

V.5  Effects  of  Normal  Stresses  on  Turbulence  Correlations 

As  noted  above  in  section  V.4  and  in  the  earlier  work  of  Simpson  et  al . 
(1977),  the  normal  stresses  turbulence  energy  production  terms  are  important  in 
separating  flows.  Simpson  et  al.  defined  a  nondimensional  factor  F  as  the 
ratio  of  total  turbulence  energy  production  to  the  shear-stress-related  turbulence 
energy  production 


-uv  aU/8y 


Figures  58  show  F-l.  Following  Collins  and  Simpson  (1976),  the  turbulence  param¬ 
eters  in  the  expression  for  F  can  be  inter-related  so  that  F  can  be  expressed  as 
a  function  of  the  rate  of  strain  ratio.  The  F  factor  can  then  be  incorporated 
into  some  of  the  turbulence  models  and  correlations  to  account  for  normal  stresses 
effects . 

Collins  and  Simpson  expressed 

(u'2  -  v'2)  =  C]7  (34) 


However,  the  present  data  available  at  a  number  of  streamwise  locations  indicate 
that  at  the  location  of  the  maximum  shearing  stress  a  better  expression  is 

(u'2-V2)^  (35) 


This  reduces  to  equation  (34)  for  flat  plate  flow  with  F  =  1.  Collins  and 
Simpson  found  C^  to  be  a  constant  equal  to  0.32  for  Klebanoff's  (1954)  zero 
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pressure  gradient  flow  and  0.28  and  0.23  for  Bradshaw's  (1967)  adverse 
pressure  gradient  flow.  In  view  of  the  definition  for  C2  in  equation  (35) 
the  separating  flow  of  Simpson  et  al .  (1977)  yields  values  for  C 2  of  0.33 
at  88  inches  where  F  =  1  and  0.44  at  103.8  and  124.3  inches.  The  distri¬ 
bution  of  C2  for  present  data  is  shown  in  Fig.  59  (b)  and  an  average  value 
of  0.375,  which  lies  within  the  experimental  uncertainty  of  26%,  was  chosen 
for  further  analysis. 

The  Reynolds  shearing  stress  can  also  be  related  to  F  and  q  by  a  modi¬ 
fication  to  Bradshaw's  correlation 

-uv  =  a2?/F4/3  (36) 

Figure  59  (c)  shows  that  this  is  a  good  fit  to  the  present  data  at  the  location 
of  the  maximum  shearing  stress  with  a2  *  0.15.  Equations  (33),  (35),  and  (36) 
can  be  combined  into  the  form 

F  =  - L -  (37) 

,  ,  l2  3U/3x 

a2  Wdy 

at  the  location  of  the  maximum  shearing  stress.  As  shown  in  Figure  59  (d), 
C2/a2  is  nearly  a  constant  within  the  experimental  uncertainty  of  +  17%  with  an 
average  value  of  2.5,  which  is  close  to  the  value  of  2.0  used  by  Collins  and 
Simpson  in  the  prediction  model  for  separating  flows.  Figure  59  (a)  shows  that 
equation  (37)  agrees  with  equation  (33)  within  the  experimental  uncertainty  of 

1 14%- 

A  two  term  binomial  expansion  of  equation  (37)  is  similar  to  Bradshaw's 
(1973)  F  factor  used  to  account  for  the  effect  of  extra  strain  rates  in  complex 
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turbulence  flows.  However,  unlike  the  case  with  Bradshaw's  factor  the  constant 
1S  derived  directly  from  the  turbulence  structure  and  is  not  just  an  empir¬ 
ical  constant  derived  from  tuning  a  prediction  method. 

As  shown  in  Figures  41  and  42,  the  mixing  length  and  eddy  viscosity  distri¬ 
butions  in  the  outer  region  decrease  in  magnitude  in  the  downstream  direction. 
This  syems  to  be  consistent  with  Gartshore's  (1967)  suggestion  of  decreased  Rey¬ 
nolds  stress  in  flows  with  an  extra  strain  rate  3V/3y,  as  in  his  own  experiments 
on  retarded  wakes.  Figures  60  (a)  and  (b)  show  these  parameters  at  the  maximum 
shearing  stress  for  each  location.  F  was  fit  to  these  data  with  the  following 
results. 


(38) 


(39) 


These  fits  were  obtained  in  the  following  manner.  The  normally  accepted 
value  of  0.08  was  used  for  ty<S  at  F  =  1.  Using  this,  an  average  value  for  it/6 
in  the  outer  region,  and  the  value  of  F  at  the  location  of  the  maximum  shearing 
stress,  the  exponenet  on  F  in  equation  (38)  was  determined  at  each  streamwise 

1  25 

location.  This  exponent  was  within  12%  of  1.25  and  the  modified  correlation  F  * 
£/6  agrees  within  the  limits  of  experimental  uncertainty  of  21%  with  the  normally 
accepted  value  of  0.08. 

For  evaluating  the  exponent  in  equation  (39),  all  values  were  taken  at  the 
location  of  the  maximum  shearing  stress;  105.3  inches  was  considered  the  location 
where  F  =  1.  Equation  (39)  agrees  with  the  data  within  the  uncertainty  of  26%. 
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V.6  Characteristic  Frequencies  from  Spectra  in  Separated  Flow 

Strickland  and  Simpson  (1975)  showed  that  the  characteristic  bursting 
frequency  could  be  determined  by  the  peak  in  the  first  moment  of  the  spectra 
nF(n)  of  the  wall  shearing  stress.  These  characteristic  frequencies  for  the 
Simpson  et  al .  (1977)  separating  flow  correlated  with  the  outer  flow  velocity 
and  length  scales,  and  6,  as  do  the  bursting  frequencies  for  the  zero- 
pressure-gradient  case.  However,  was  between  11.7  and  8.35  for  that 

separating  flow,  whereas  values  of  about  5  are  reported  for  the  zero-pressure- 
gradient  case.  The  basic  conclusion  of  these  earlier  results  is  that  the 
characteristic  frequency  of  the  most  energetic  turbulent  fluctuations  scale  on 
the  large-scale  structure  of  the  shear  flow. 

In  the  earlier  work  of  Simpson  et  al.  (1977)  no  spectral  measurements  in 
the  separated  flow  were  made.  In  the  present  flow  spectral  data  for  u  were 
obtained  from  the  laser  anemometer  velocity  signals.  Since  the  LDA  signal  data 
rate  was  under  400  signals  per  second  and  signal  dropout  was  present,  the  spectra 
are  only  reliable  under  100Hz.  The  first  moment  of  each  spectral  distribution 
nF(n)  was  obtained  and  the  frequency  of  the  peak  was  selected  as  the  characteristic 
frequency  nb.  In  many  cases  the  nF(n)  peak  was  constant  over  a  frequency  range, 
which  is  represented  in  Figures  61  as  a  line  over  the  range  of  U^/Sn^  values  for 
a  spectrum  at  a  given  y/6. 

Figure  61a  shows  that  upstream  of  intermittent  separation  U^/Snh  is  essenti¬ 
ally  constant  throughout  the  inner  flow  region  with  a  value  of  about  10  +  3.  At 
successive  downstream  locations  the  range  of  U^/Sn^  for  a  given  nF(n)  peak  be¬ 
comes  progressively  larger  near  the  wall  as  shown  in  Figure  61  (b-f).  In  most 
cases  a  single  frequency  characterises  the  nF(n)  peak  in  the  outer  region.  U^/finh 
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Figure  61(b). 
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is  about  10+3  at  the  lower  end  of  the  U^/Sn^  bands.  The  upper  values 
of  uySn^  are  about  40  or  so  in  the  inner  region. 

These  results  indicate  that  the  characteristic  frequency  of  the  outer 
region  correlate  with  and  6  along  the  flow,  with  an  approximately  constant 
value  of  U  /6n.  of  about  10  +  5.  This  is  consistent  with  the  earlier  work  of 
Simpson  et  al.  (1977).  Nearer  the  wall  the  frequency  range  of  the  energetic 
turbulent  motions  descends  to  frequencies  one-fourth  as  large. 

For  attached  boundary  layers  the  spectra  for  the  near  wall  flow  have  a 
range  of  frequencies  over  which  the  peak  of  each  nF(n)  distribution  is  constant 
(Rotta,  1962).  This  is  a  consequence  of  the  logarithmic  law-of-the-wall  velocity 
profile.  For  a  separated  flow  the  law-of-the-wall  is  not  valid,  so  a  different 
explanation  of  the  nF(n)  distribution  near  the  wall  is  needed.  The  upper  fre¬ 
quency  end  of  the  nF(n)  peak  is  at  approximately  the  same  frequency  as  the  outer 
region  peak  frequency.  Note  from  Figs.  61  (b-f)  that  the  wide  frequency  spectral 
peaks  seem  to  occur  at  locations  near  the  wall  where  yu  <  1. 

One  simple  speculation  is  that  the  celerity  or  speed  of  the  eddies  in  the 
backflow  region  is  much  lower  than  that  in  the  outer  region.  Fig.  17  of  Simpson 
et  al.  (1977)  supports  this  idea.  Thus,  as  large  scale  structures  pass  through 
the  outer  flow  at  a  frequency  of  about  U^/106,  these  same  structures  move  at  a 
much  lower  average  celerity  in  the  backflow  region,  producing  a  much  lower  fre¬ 
quency  spectrum. 

VI.  CONCLUSIONS  -  The  Nature  of  a  Separating  Turbulent  Boundary  Layer 

These  experiments  confirm  the  conclusions  of  Simpson  et  al.  (1977)  re¬ 
garding  a  separating  airfoil  type  turbulent  boundary  layer.  The  mean  flow  up¬ 
stream  of  the  beginning  of  intermittent  separation  obeys  the  law-of-the-wall  and 
the  Perry  and  Schofield  (1973)  velocity  profile  correlation  for  the  outer  region. 
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Sandborn's  correlations  for  the  locations  of  intermittent  separation  (y  = 

0.8)  and  fully-developed  separation  hold.  Pressure  gradient  relaxation  begins 
upstream  of  intermittent  separation  near  the  wall  jet  control  in  this  flow  and 
continues  until  the  location  of  fully-developed  separation.  The  upstream- 
downstream  intermittency  y  ,  u',  and  v',  and  -u7  profiles  each  approach  similar¬ 
ity  profiles  downstream  of  separation.  The  frequency  of  passage  of  the  outer 
region  large  scale  eddies  n^  scales  on  the  free-stream  velocity  Uto  and  the  boundary 
layer  thickness  6.  Normal  stresses  effects  contribute  significantly  to  the  momen¬ 
tum  and  turbulence  energy  equations. 

Much  new  information  about  the  senarated  region  has  been  gathered  and  leads 
to  significant  conclusions  about  the  nature  of  the  separated  flow.  For  reference 
the  most  important  results  are  summarized  below. 

1.  The  backflow  mean  velocity  profile  scales  on  the  maximum  negative  mean 

+  4. 

velocity  and  its  distance  from  the  wall  N.  A  U  vs.  y  law-of-the- 
wall  velocity  profile  is  not  consistent  with  this  correlation  since 
both  and  N  increase  with  streamwise  distance,  while  the  law-of-the- 
wall  length  scale  v/UT  varies  inversely  with  the  velocity  scale  U . 

It  does  not  appear  possible  to  describe  the  separated  flow  mean  velocity 
profiles  by  a  universal  "backflow  function"  that  is  added  to  a  universal 
"wake  function". 

2.  High  turbulence  levels  exist  in  the  backflow,  u'  and  v'  are  of  the 
same  order  as  Ju|  .  Since  the  free-stream  velocity  in  the  separated 
region  is  rather  steady,  this  means  that  the  near  wall  fluctuations 
are  not  mainly  due  to  a  flapping  of  the  entire  shear  layer,  but  are  due 
to  turbulence  within  the  separated  shear  layer. 
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3.  Low  levels  of  Reynolds  shearing  stress  occur  in  the  backflow. 

-uv/uV  and  -uv/(u'  +  v'  )  correlations  are  low  in  the  backflow, 

but  are  comparable  with  those  for  mixing  layers  in  the  outer  region. 

4.  Mixing  length  and  eddy  viscosity  models  are  adequate  upstream  of 
separation  and  in  the  outer  region,  but  are  physically  meaningless 
in  the  backflow.  Normal  stresses  effects  appear  to  account  for  the 
lower  mixing  length  and  eddy  viscosity  values  observed  in  the  outer 
region  of  the  separated  flow. 

5.  In  the  separated  flow  between  the  wall  and  the  locations  of  the 
minimum  mean  velocity,  the  skewness  factor  for  u,  Su,  is  negative. 
Between  this  point  and  the  locations  of  the  maximum  shearing  stress, 

Su  is  positive.  The  flatness  factor  reaches  a  local  maximum  of 
about  4  at  the  minimum  mean  velocity  location.  Sy  has  a  profile 
shape  and  magnitudes  that  are  approximately  the  mirror  image  or 
negative  of  Su. 

6.  Negligible  turbulence  energy  production  occurs  in  the  backflow. 

Normal  and  shear  stresses  production  in  the  outer  region  supply 
turbulence  energy  to  the  backflow  by  turbulent  diffusion.  These 
results  are  consistent  with  the  absence  of  a  logarithmic  mean  velocity 
profile  in  the  backflow,  since  classical  turbulence  energy  production 
arguments  indicate  that  the  rate  of  production  must  equal  the  rate  of 
dissipation  in  such  a  region. 

These  turbulence  energy  results  lead  to  the  conclusion  that  the  backflow 
is  controlled  by  the  large-scale  outer  region  flow.  Movies  of  laser-illuminated 
smoke  also  have  clearly  revealed  that  the  large  eddy  structure  supplies  most  of 


near  wall  backflow.  The  small  mean  backflow  does  not  come  from  far  down¬ 
stream  as  suggested  in  Figure  62(a),  but  appears  to  be  supplied  intermittently 
by  large-scale  structures  as  they  pass  through  the  separated  flow  as  suggested 
by  Figure  62(b). 

A  simple  qualitative  experiment  was  performed  to  determine  qualitatively 
the  influence  of  the  downstream  near  wall  conditions  on  the  separation  behavior. 

A  deflection  plate  was  located  at  the  end  of  the  second  section  (200  inches) 
as  shown  in  Figure  63.  For  heights  of  this  deflection  plate  up  to  7  inches, 
no  appreciable  change  in  the  separation  zone  location  (122-140  inches)  and 
behavior  were  noted.  This  result  also  seems  to  support  the  flow  model  suggested 
in  Figure  62(b)  where  the  backflow  is  supplied  locally  by  outer  region  large- 
scale  structures.  Only  after  the  deflection  plate  was  high  enough  to  begin  to 
change  the  free-stream  pressure  gradient  did  the  location  of  the  separation 
zone  change. 

Of  course,  this  mechanism  for  supplying  the  backflow  may  be  dominant  only 
when  the  thickness  of  the  backflow  region  is  small  as  compared  with  the  turbulent 
shear  layer  thickness,  as  in  the  present  case.  Experiments  (Fox  and  Kline,  1962) 
on  separation  in  wide-angle  diffusers  indicate  that  the  mean  backflow  can  come 
from  far  downstream  when  the  thickness  of  the  backflow  region  is  comparable  to 
the  thickness  of  the  forward  flow. 

Downstream  of  fully-developed  separation  in  these  experiments,  the  mean 
backflow  region  appears  to  be  divided  into  three  layers:  a  viscous  layer  nearest 
the  wall  that  is  dominated  by  the  turbulent  flow  unsteadiness  but  with  little 
Reynolds  shearing  stress  effects;  a  rather  flat  intermediate  layer  that  seems  to 
act  as  an  overlap  region  between  the  viscous  wall  and  outer  regions;  and  the 
outer  backflow  region  that  is  really  part  of  the  large-scaled  outer  region  flow. 


zero  average  wall  shear  stress 

Turbulent  boundary  layer  Separated  flow  region 


Figure  62  (a).  Traditional  view  of  turbulent  boundary  layer  separation  with  the 

MEAN  BACKFLOW  COMING  FROM  FAR  DOWNSTREAM. 


Figure  62  (b),  A  possible  flow  model  with  the  coherent  structures  supplying 

THE  SMALL  MEAN  BACKFLOW. 


173 


The  Reynolds  shearing  stresses  in  this  region  must  be  modeled  by 
relating  them  to  the  turbulence  structure  and  not  to  local  mean  velocity 
gradients.  The  mean  velocity  profiles  in  the  backflow  are  a  result  of 
time-averaging  the  large  turbulent  fluctuations  and  are  not  related  to  the 
cause  of  the  turbulence.  In  contrast,  in  flows  for  which  the  eddy  viscosity 
and  mixing  length  models  appear  to  be  useful,  the  instantaneous  velocity 
gradients  are  not  extremely  different  from  the  local  mean  velocity  gradient 
and  significant  local  turbulence  energy  production  occurs,  i.e.,  the  Reynolds 
shearing  stresses  is  physically  related  to  the  mean  velocity  gradient. 

VII.  FUTURE  WORK 

Currently  measurements  of  w'  and  ^  are  being  made  in  the  separated 
flow  to  completely  document  this  flow.  During  the  1980-81  period  a  scanning 
laser  anemometer  system  will  be  developed  to  obtain  almost  instantaneous 
velocity  profiles.  These  instantaneous  profiles  should  prove  useful  in 
relating  the  instantaneous  backflow  to  the  outer  region  flow. 
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APPENDIX 


TABULATION  OF  LASER, 

CROSS  HOT-WIRE  AND  SINGLE 
HOT-WIRE  ANEMOMETER  DATA 

(Note  that  only  first  three  digits  in 
each  number  are  valid.) 
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